Near-infrared spectroscopy of five Blue Compact Dwarf galaxies: 
II Zw 40, Mrk 71, Mrk 930, Mrk 996 and SBS 0335-052Ei 



Yuri I. Izotov 

Main Astronomical Observatory, Ukrainian National Academy of Sciences, 27 Zabolotnoho 

str., Kyiv 03680, Ukraine 

izotovQmao . kiev . ua 

and 

Trinh X. Thuan 

Astronomy Department, University of Virginia, P.O. Box 400325, Charlottesville, VA 

22904-4325 

txtSvirginia . edu 
ABSTRACT 

We present near-infrared (NIR) spectroscopic observations of five blue com- 
pact dwarf (BCD) galaxies, II Zw 40, Mrk 71, Mrk 930, Mrk 996 and SBS 0335- 
052E. The NIR spectra which cover the 0.90 fim - 2.40 fim wavelength range, 
show hydrogen, molecular hydrogen, helium, sulfur and iron emission lines. The 
NIR data for all BCDs have been supplemented by optical spectra. We found 
the extinction coefficient in all BCDs to be very similar in both the optical and 
NIR ranges. The NIR hydrogen emission lines do not reveal more star forma- 
tion than seen in the optical. The same conclusion is reached from Spitzer data 
concerning the MIR emission lines. This implies that emission-line spectra of 
low-metallicity BCDs in the ~ 0.36 - 25 /im wavelength range are emitted by 
relatively transparent ionized gas. The large extinction derived from the MIR 
continuum emission in some BCDs implies that the latter arises not from the vis- 
ible H II regions themselves, but from locations outside these H ll regions. The 
H2 emission line fluxes can be accounted for by fluorescence. CLOUDY stellar 
photoinization models of all BCDs reproduce well the fluxes of most of the ob- 
served optical and NIR emission lines, except in Mrk 930 where shock ionization 
is needed to account for the [Fe ll] emission lines. However, some contribution 
of shock ionization at the level of < 10% that of stellar ionization is required to 
reproduce the observed fluxes of high ionization species, such as He ll AO. 469 /im 
in the optical range and [O iv] A25.89 /im in the MIR range. 
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Subject headings: galaxies: abundances — galaxies: irregular — galaxies: ISM 
— H II regions — galaxies: individual (II Zw 40, Mrk 71, Mrk 930, Mrk 996, 
SBS 0335-052E) — infrared: ISM 



INTRODUCTION 



Blue compact dwarf (BC D) galaxies are the least chemically evolved star-forming galax- 
ies known in the universe (see iThuan II2008I . for a review). They have an oxygen abundance 
in the ra nge 12 + log 0/H = 7.1 - 8.3, i.e. 1/40 -1/3 that of the Sun if the solar abun- 
dance of lAsplund et al.l ( 120091 ). 12 -|- log 0/H = 8.7, is adopted. They thus constitute 



excellent nearby laboratories for studying star formation and its interaction with the in- 
terstellar medium (ISM) in a very metal-deficient environment, intermediate between the 
solar-metallicity environment of galaxies such as the Milky Way, and the pristine environ- 
ment of primordial galaxies. While BCDs have been extensively studied spectroscopically 
in the optical range, not much work has been done in the near-infrared (NIR) 0.9 - 2.4 fim 
range. We have thus started a long-term observational program to study the NIR spectro- 
scopic properties of a sample of BCDs chosen to cover a large range of metallicities, with 12 
-|- log 0/H ranging from 7.1 to 8.3. The first result s of the program, co ncerning the BCD 
Mrk 59 (12 + log 0/H = 8.0), have been discussed by llzotov et al. I (120091 ). Here, we present 
new J, H and K spectroscopic observations of five more BCDs, with 12 -|- log 0/H ranging 
from 7.29 to 8.10. 



In this new 
been studied by 



ist, only two B CDs h ave had previous NIR obse rv ations. II Zw 40 ha s 



Vanzi et al. I (119961 ). ICoziol et al. 



has been observed by 



Vanzi et al. 



,Tov fc Lester I (119881) . iMoorwood fc Oliva I (119881) . iMouri et al. I (119901 ) 



20011) and IVanzi et al. I (120081), while SBS 0335-052E 



J2OO0I ). However, all these observations covered a more 
restricted wavelength range, including only either the H band and/or the K band. The 
other three BCDs, Mrk 71, Mrk 930 and Mrk 996, have not been observed spectroscopi- 
cally in the NIR before. The general characteristics of the five objects are given in Table 
dJ With advances in NIR detector sensitivities, and with the commisionning in 2008 of the 
TripleSpec spectrograph at the Apache Point Observatory (APO) 3.5-meter telescop^, we 
can obtain NIR spectra with a simultaneous coverage of the JHK bands and considerably 



^Tables [3] -[7] are presented only in the electronic edition of the Journal. 

^The Apache Point Observatory 3.5-nieter telescope is owned and operated by the Astrophysical Research 
Consortium. 
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more signal-to-noise. These improved observations will allow us to study in more detail the 
physical conditions in BCDs, and check for their dependence on metallicity. We will be able 
to compare the extinction in the optical and NIR and search for hidden star formation. We 
will also be able to study the excitation mechanisms of line emission of ionized species in the 
H II region, and search for molecular hydrogen emission. 

In addition to the NIR observations, we have also obtained optical spectra of the same 
objects, except for Mrk 930, in the wavelength range 0.36 - 0.97 ^m. These were taken 
so to have an overlap with the NIR spectra in the 0.90 - 0.97 /xm range, which allows us 
to match the two types of spectra and reduce the uncertainties due to different apertures 
and/or varying weather conditions. 

We describe the 3.5m APO observations for the five new BCDs in Section [21 In Section 
[31 we discuss their extinction, both in the optical and NIR ranges, and in the MIR range 
using Spitzer data, the excitation mechanisms for molecular hydrogen and iron emission in 
the NIR range, and CLOUDY modeling of the H ll region. We summarize our findings in 
Section [H 



2. OBSERVATIONS 

2.1. Near-infrared observations 

NIR spectra of the BCDs II Zw 40, Mrk 71, Mrk 930, Mrk 996 and SBS 0335-052E were 
obtained with the 3.5 m APO telescope in co njunction with the T ripleSpec spectrograph on 



the night of 2009 November 25. TripleSpec (jWilson et al. II2004I ) is a cross-dispersed NIR 



spectrograph that provides simultaneous continuous wavelength coverage from 0.90 to 2.46 
/im in five spectral orders during a single exposure. A l'.'lx43" slit was used, resulting in a 
resolving power of 3500. 

Three AOV standard stars, HD 16286, HD 23504, HD 212170 were observed for flux 
calibration and correction for telluric absorption. Spectra of Ar comparison arcs were ob- 
tained on the same date for wavelength calibration. Since all observed targets are smaller 
than the length of the slit, the nod-on-slit technique was used to acquire the sky spectrum. 
Objects were observed by nodding between two positions A and B along the slit, following 
the ABBA sequence, and with an integration time of 200s or 300s at each position. The 
journal of the NIR observations is given in Table [2l We show there the observation date, 
the exposure time, the extraction aperture of the one-dimensional spectrum and the average 
airmass. 
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We carried out the reduction of the data according to the following procedures. The 
two-dimensional spectra were first cleaned for cosmic ray hits using the IRAFi routine CR- 
MEDIAN. Then all A and B frames were separately coadded and the resulting B frame 
was subtracted from the resulting A frame. Finally, the (negative) spectrum at position B 
was adjusted to the (positive) spectrum at position A and subtracted from it. The same 
reduction scheme was applied to the standard stars. 

We then use the IRAF routines IDENTIFY, REIDENTIFY, FITCOORD, TRANS- 
FORM to perform wavelength calibration and correction for distortion and tilt for each 
frame. For all galaxies, a one-dimensional spectrum was extracted from the two-dimensional 
frames using the APALL IRAF routine, within a I'.'lxG" rectangular aperture so as to in- 
clude the brightest star-forming regions. The redshift of our objects range from 0.00036 
for Mrk 71 to 0.0183 for Mrk 930, so that 1" corresponds to 7 pc for the nearest object 
and 350 pc for the furthest object, adopting a Hubble constant of 75 km s~^ Mpc~^. For 
the nearer objects, II Zw 40, Mrk 71 and Mrk 996, the slit contains a single stellar cluster 
which dominates in the line emission. Thus, in II Zw 40 the slit is centered on the compact 
su per-star cluster (SSC) A (IVanzi et al. II2008I ). in Mrk 71 on the brightest compact cl uster 
A ( iThuan fc Izotov II2005I ). and in Mrk 996 on the central cluster (IThuan et al. Ill996l ). All 
these stellar clusters are unresolved on HST/WFPC2 images. For the furthest galaxies, SBS 
0335-052E and Mrk 930, the slit is also centered on the brightest cluster, but contains sev- 
eral other stellar clus ters. The HST images of Mrk 930 and SBS 0335-052E are shown in 
Malkan et alTI Jl998l ) and iThuan et al. I f ll997h . respectively. For SBS 0335-052E, the slit 
is oriented at a position angle of -30° and includes not only SSCs 1 and 2, which dominate 
the line emission, but also SSCs 3-6. The spectrum of of Mrk 930 (its southern knot) was 
obtained at a position angle of +238°. Because of the compact structure of SBS 0335-052E 
and Mrk 930, we have not attempted to extract spectra of individual clusters and have used 
the integrated spectra. For these two furthest objects, the larger linear sizes of the slit also 
mean that a larger area of the weakly ionized and neutral ISM around the stellar clusters 
is included in the slit, a fact which will be relevant in our modelization of the line emission 
of the low- ionization and neutral species (Section 13. 3p and in our discussion of the MIR 
extinction (Section l3.6p . 

Flux calibration and correction for telluric absorption were performed by first multiply- 
ing the one-dimensional spectrum of each galaxy by the synthetic absolute spectral distri- 
bution of the brightest standard star HD 23504, smoothed to the same spectral resolution. 



•^IRAF is the Image Reduction and Analysis Facility distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under 
cooperative agreement with the National Science Foundation (NSF). 
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and then by dividing the result by the observed one-dimensional spectrum of the same star. 
Since there does not exist any published absolute spectral energy distribution of standard 
stars, we have derived the one of HD 23504 by scaling the synthetic absolute spectral energy 
distribution of the star Vega [a Lyrae), of similar spectral type AOV, to the brightness of 
HD 23504. 

The resulting one- dimensional flux-calibrated and redshift-corrected NIR spectra of all 
five BCDs are shown in Figs. [U-El 



2.2. Optical observations 

To have a more complete physical picture of the extinction and of the ionization mecha- 
nisms in our five BCDs, we have supplemented the NIR observations with optical ones. The 
optical observations of all galaxies, with the exception of Mrk 930, were obtained with the 3.5 
m APO telescope in the course of several nights during the period N ovember 2008 - Februar y 



2010. For Mrk 930, we have used the optical spectrum obtained by llzotov &: Thuan I ( 119981 ). 



The 3.5 m APO observations were made with the Dual Imaging Spectrograph (DIS) 
in both the blue and red wavelength ranges. In the blue range, we use the B400 grating 
with a linear dispersion of 1.83x10"^ /im/pix and a central wavelength of 0.44 fim, while 
in the red range we use the R300 grating with a linear dispersion 2.31x10"^ /xm/pix and a 
central wavelength of 0.75 /xm. The above instrumental set-up gave a spatial scale along the 
slit of 0'.'4 pixel"^, a spectral range ~ 0.36 - 0.97 /im and a spectral resolution of 7x10""^ 
fim (FWHM). Several Kitt Peak IRS spectroscopic standard stars were observed for flux 
calibration. Spectra of He-Ne-Ar comparison arcs were obtained at the beginning or the end 
of each night for wavelength calibration. The journal of the optical observations is given in 
Table [21 As for Mrk 930, the acquisition of its optical spectrum in the wavelength range ~ 



0.36 - 0.74 /xm is described in llzotov &: Thuan 1 (119981 ). 



The two-dimensional optical spectra were bias subtracted and fiat-field corrected using 
IRAF. We then use the IRAF software routines IDENTIFY, REIDENTIFY, FITCOORD, 
TRANSFORM to perform wavelength calibration and correction for distortion and tilt for 
each frame. One-dimensional spectra were then extracted from each frame using the APALL 
routine. The extraction apertures for the one-dimensional optical spectra are shown in Table 
121 The areas of these optical extraction apertures are generally similar to those of the NIR 
extraction apertures. An exception is Mrk 930, for which the area of the optical extraction 
aperture is ~ 2 times that of the NIR extraction aperture. The sensitivity curve was obtained 
by fitting with a high-order polynomial the observed spectral energy distribution of standard 
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stars. Then sensitivity curves obtained from observations of different stars during the same 
night were averaged. The flux deviation from curve to curve is less than 1-2%. 

We have used the optical spectra to derive the physical c onditions and oxygen abun- 



dances of all five BCDs, following the procedures described by llzotov et al. I (l2006al ). The 
derived oxygen abundances are shown in Table [H 

The NIR and optical spectra of all BCDs, except for Mrk 930, overlap in the 0.90 - 

0. 97 /xm wavelength range, allowing us to match the fluxes of strong emission lines that are 
common in the two types of spectra. Specifically, we have adjusted the optical and NIR 
spectra so that the flux of the strongest line in the overlapping region, [S ill] 0.953 /im, is the 
same in both spectra. We have thus multiplied the NIR spectra of II Zw 40, Mrk 71, Mrk 
996 and SBS 0335-052E by factors of 1.27, 1.17, 1.10 and 1.09, respectively. The multiplying 
factors are close to unity, as expected from the similarities of the areas of the optical and 
NIR extraction apertures. The deviations from unity can be caused by the slighly differing 
extraction apertures, non-perfect pointing and varying weather conditions (transparency, 
seeing). They can also be caused by uncertainties in the flux calibration of the NIR spectra 
introduced by uncertainties in the absolute spectral energy distribution (SED) of the NIR 
standard star HD 23504. As discussed in Sect. 12. H this SED is not known and was obtained 
by scaling the synthethic NIR spectrum of Vega. 

Additional uncertainties are introduced by the ionization structure of the H ll regions, 

1. e. differences in the spatial distribution of various ions. In particular, hydrogen line emission 
is more concentrated to the center of the star-forming region than the [S ill] emission. We 
find that the FWHM of the Ha emission along the spatial axis of our objects is ~20% 
smaller than that of the [S ill] A0.953/im emission. The ionization structure uncertainties 
are difficult to estimate as the observations were obtained with different slit widths, 1'.'5 for 
optical spectra and I'.'l for NIR ones. We have attempted to minimize these uncertainties 
by extracting one- dimensional spectra within apertures with similar areas. 

For Mrk 930, we have scaled the NIR spectrum by a factor of 1.73 so that the Bvy/HP 
flux ratio is equal to the theoretical recombination value, after correcting the emission line 
fluxes for interstellar extinction and underlying stellar absorption derived from the optical 
spectrum. This scaling factor agrees nicely with the ratio of the areas of the apertures used 
for the NIR and optical observations of Mrk 930. The combined redshift-corrected optical 
and NIR spectra for all BCDs are shown in Figs. [6]-[T0l Gaps in the NIR spectra correspond 
to regions of strong absorption by telluric lines. While the scaling of the NIR spectrum of 
Mrk 930 is more uncertain, it seems justifled by the fact that NIR SED extends smoothly 
the optical SED, without any abrupt discontinuity (Fig. [8]). 
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RESULTS AND DISCUSSION 



Emission-line fluxes in both the optical and NIR ranges were measured by using the 
SPLOT routine in IRAF. The errors of the line fluxes were calculated from the photon 
statistics in the non-flux-calibrated spectra. They do not include systematic errors such 
as deviations of the standard star SED from the Vega SED, different apertures used in 
the optical and NIR observations and varying weather conditions. Those are difficult to 
estimate as discussed in Sect. 12.21 The observed fluxes F{X) of emission lines derived from 
the adjusted optical and NIR spectra are shown in the second columns of Tables [3] - [71 
They are normalized to the observed H/3 flux, F(H/3), shown at the end of each table, and 
multiplied by a factor of 100. 

How do our ob servations coni p are w ith p revious ones o f the same objects? II Zw 
40 was observed by IJov fc Lester I (119881) and ICoziol et al. I (120011) i n the K band, and 



by iMoorwood &: Oliva I ( 



K bands. 



Jov &: Lester 



9881) ■ 



(Il988r ) (3'.'6 aperture). 



Vanzi et al. I (119961) and IVanzi et al. (120081 ) in the H and 



Vanzi et al. (jl996 ) (2'.'4 aperture) and 



Coziol et al. I ( 20011 ) (2-3" aperture) obtained a flux of the Br7 emission line F(Br7) in 
the range (3.4-3.9) xlO~^^ erg s^^ cm^^, slight ly larger than our va l ue of 3.1x10"^^ erg s~^ 
cm^^ (l'.'lx6" aperture). On the other hand, iMoorwood fc Oliva 



6.7x10 erg s ^ cm ^ in a 6"x6" aperture, while Vanzi et al. ( 20081 ) derived F(Br7) 



1119881) found F(Br7) 



1.0x10 erg s ^ cm ^ with in a circular apertur e of 1" in radius. SBS 0335-052E has been 

1 ^ —1 
erg s 



observed in the NIR only by IVanzi et al. I ( 120001 ) who derived F(Br7) = 1.5x10 
cm~^ in an aperture of l"xl'.'5, lower than our value of 2.4x10"^^ erg s~^ cm"^. All ob- 
servations of II Zw 40 and SBS 0335-052E appear to be consistent in the sense that the 
Br7 fluxes obtained with apertures having similar areas are close to each other while those 
obtained with smaller apertures are lower. 

The presence of many emission lines in the optical and NIR spectra (Figs. HI - HOl and 
Tables 131-171) allows us to study extinction, H2 and [Fe ll] emission and excitation mechanisms 
in all five galaxies. We discuss these issues in the following sections. 



3.1. Optical and NIR extinction and hidden star formation 



^In previous optical- NIR studies of BCDs, except for our own work on Mrk 59 ([izotov et al. 

20091 ). JHK spectra have been observed separately, and there has been no wavelength over- 



lap between the optical and NIR spectra. This introduces uncertainties, due to the use of 
generally different apertures in the optical and NIR observations, and due to the adjusting of 
the continuum levels of the NIR spectra obtained separately in different orders. Because our 
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NIR spectra have been obtained simultaneously over the whole JHK wavelength range and 
because it is possible to directly match the optical and NIR spectra using common emission 
line fluxes (except for Mrk 930), these adjusting uncertainties are mostly eliminated. This 
allows us to directly compare the optical and NIR extinctions and search for any hidden 
extinction at longer wavelengths. 

In the second column of Tables [3] - [TJ we list the observed fluxes F{X) of the emission 
lines in both the optical and NIR spectra. In the third column of the same Tables are shown 
the fluxes /(A), corrected for underlying stellar absorption (for hydrogen lines) and extinction 
(for all lines). The extinction coefficient C(H/3) and equivalent width EW(abs) of hydrogen 
lines derived from the decrement of hydrogen Balmer lines in the optical spectrum are listed 
in the footnotes to Tables E]- El They are used to correct the emission line flux es in both the 



optical and NIR ranges, adopting the extinction curve of ICardelli et al. I ( 119891 ) and RiV) = 
A{y) / E{B — V) = 3.2, where A{y) and E{B — V) are respectively the V band extinction 
and the B — V reddening. The value R{y) =3.2 is typical for the diffuse interstellar medium 
in the Galaxy. With this value of R{V), then C(H/3) = 1.47x^(fi -V) = 0.46xA(\/). 

However, the dust composition in the low-metallicity interstellar medium in our BCDs 
can be different from tha t in th e Galaxy, resulting in a different extinction curve. Indeed, 



Rocca-Volmerange et al. I ( 119811 ) found that the UV extinction curve in the Small Magellanic 
Cloud (SMC) is signiflcantly steeper than that in the Galaxy, indicating the presence of a 
larger fraction of small dust grains in the SMC. To reproduce the SMC extinction curve in 
the UV range, the value of RiV) should be smaller than 3.2. However, differences in the 
extinction curves with different R{y)s are much smaller in the optical and NIR ranges. For 
example, if we adopt RiV) = 2.5, corresponding to C(H/3) = 1.17xE{B -V) = 0A7xA{V), 
then C(H/3) is decreased by ~ 10% and the flux of the Br7 emission line relative to H/3 is 
increased only by ~ 2 - 3%, much smaller than the errors given in Tables [31 - 13 for that line. 
Therefore, in the following, we will use the extinction-corrected emission-line fluxes obtained 
with R{V) = 3.2. 

The corrected fluxes of hydrogen lines in the third column of the same Tables are 
derived, assuming that the difference between the observed and theoretical hydrogen Balmer 
decrements is solely due to extinction and underlying stellar absorption. However, collisional 
excitation of hydrogen lines may play a role, resulting in a steepening of the observed Balmer 



decrement. Following iLuridiana I ( 120091 ). we assume that the fraction of the Ha emission due 
to collisional excitation is 5%, while that for H/3, H7 and H/3 emission is 3%. These values 
should be considered as upper limits for c ollisional excitation. A possible exception is Mrk 



996 with a very high-density H ll region (iThuan et al. 1120081 ). where collisional excitation 



could be higher. Subtracting the collisional contribution from the Ha, H/3, H7 and B.6 
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emission lines and correcting all lines for extinction and hydrogen lines for underlying stellar 
absorption, we obtain the emission line fluxes listed in the fourth column of Tables |3] - [71 

In the second to sixth columns of Table [HI we show the extinction-corrected fluxes /(A) 
of the strongest hydrogen lines in the optical and NIR spectra for all observed galaxies. For 
comparison, the seventh column of Tab le [8] shows the theoretical recombination flux ratios, 
calculated by [Hummer fc Storey I (1l987l ) for case B, with an electron temperature Tg = 15000 
K and an electron number density A'^e = 100 cm~^. It can be seen that there is a very good 
agreement between the corrected and theoretical recombination values of the optical and 
NIR line fluxes. We have shown in Table [9] the extinction coefficients C(H/3) derived from 
the ratio of the flux of each individual hydrogen line to the H/3 flux. They were derived from 
the equation 



/(A) F(A) EW{X) + EW{abs) 



EW{}1(3) 



/(H/3) F(H/3) 



EW{X) 



EW{RI3) + EW{abs) 



^ ;^q/(A)C(H/3) 



where /(A) is the value of the normalized reddening curve from [Cardelli et al. I ( 1l989[ ) at 
the wavelength A, adopting A{V)/E{B — V) = 3.2, F and / are respectively the observed 
and theoretical fluxes, EW(A) and EW(H/9) are equivalent widths of emission lines, and 
EW(abs) is the equivalent width of the absorption hydrogen lines. We adopt EW (abs) = 
2.7A, corresponding to a starburst age of < 4 Myr ( Gonzalez Delgado et al.lll999 ). 



Inspection of Table [9] shows that, for each object, the extinction coefficients C(H/3) 
derived from different hydrogen lines agree within the errors and do not show any systematic 
trend. This implies that the same extinction coefficient C(H/3) [or extinction ^(^)] holds 
over the whole 0.36 - 2.40 //m wavelength range. The fact that A{V) does not increase at 
longer wavelengths implies that there is no hidden star formation in the NIR as compared to 



in the optical range. Th is appears to be a general result for BCDs (e.g. [Vanzi et al. I [2000 



200a [izotov et al. I[2009[ ). This is true even in the extreme case of the BCD SBS 0335-052E 



where observations of the far-infrared (FIR) continuum imply a very large extinction [v4(y) 



(Thuan et al. 


1999; 


Houck et al. 


2004) 



3.2. H2 emission 

Molecular hydrogen lines do not originate in the H ll region, but in neutral molecular 
clouds. In the NIR, they are excited by two mechanisms. The first one is the thermal mech- 
anism consisting of collisions between neutral species (e.g., H, H2), resulting from large-scale 
shocks driven by powerful stellar winds, supernova remnants or molecular cloud collisions. 
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The second one is the fluorescent mechanism due to absorption of ultraviolet photons. It is 
known that these two mechanisms excite mostly different roto- vibrational levels of H2. By 
cor aparing the observed line ratio s with those predicted by models such as those calculated 
by iBlack &: van Dishoeck I ( Il987l ). it is possible to discriminate between the two processes. 
In particular, line emission from the vibrational level v=2 and higher vibrational levels are 
virtually absent in coUisionally excited spectra, while they are relatively strong in fluorescent 
spectra. 

H2 emission lines are detected in the NIR spectra of all our objects (see Figs. 1-5). 
Only one H2 emission line, 2.122 /im 1-0 S(l), is definitely detected in the spectrum of 
SBS 0335-052E. Another line, 2.2 48 /xm 2-1 S ( l), is marginally detected, although it is 
clearly present in the spectrum of Vanzi et al. (j2000 ). Three, five, thirteen and fourteen 
H2 lines are respectively detected in the NIR spectra of Mrk 930, Mrk 996, II Zw 40 and 
Mrk 71. The fluxes of the H2 lines relative to that of the strongest 2.122 /im 1-0 S(l) line 
are shown in Table [TOl In the la st two columns are shown the theoretical ratios calculated 
by iBlack fc van Dishoeck I ( 119871 ) in the cases of fluorescent and collisional excitation. It is 
seen that the observed line flux ratios are in agreement with those predicted for fluorescent 
excitation, despite the rather high errors in the flux ratios. Ou r flnding is i n agre ement 



with the conc lusions of I Vanzi et al. I fcOOOf ) for SBS 0335-052E, of I Vanzi et al. I mm for II 



Zw 40 and of llzotov et al. I ( 120091 ) for Mrk 59, where it was also found that the fluorescence 
process is the main excitation mechanism of H2 lines. This conclusion appears to be common 
for BCDs with high-excitation spectra. 



3.3. CLOUDY stellar photoionization modeling of the H ll region 



We next examine the excitation mechanisms of the NIR emission lines arising in the 
H II regions of our BCDs. For this purpose, we have constructed stellar photoion i zation 
mo dels for the H ll region s using the CLOUDY code (version 07.02.01) of iFerland I (Il996[ ) 



and iFerland et al. I ( 119981 ). We list in Table ITT] the input parameters of these models. The 
flrst row gives (5(H), the logarithm of the number of ionizing photons per second, calculated 
from the H/3 luminosity with the distance given in Tabled! For the determination of the H/3 
luminosity, we have taken the H/3 fluxes from large aperture spectroscopic and photometric 
observations. We t hus adopt the H/3 fluxes to be 3.5xl0 ~^^ erg s~^ cm~^ for II Zw 40 
(Ihagos et aini2007[ l. 2.1x10"^^ erg s"^ cm"^ for Mrk 71 JMoustakas k Kennicutt I [20061 ) 



and 2.5x10 erg s ^ cm ^ for Mrk 930 ( Moustakas fc Kennicutt I 20061 ). These fluxes were 



then corrected for extinction, as derived from the spectroscopic observations. We also adopt 



the extinction-corrected H/3 fluxes of 6.6x10 erg s cm for Mrk 996 ( iThuan et al. 
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20081) and of 1.8xl0~^=^ erg s"^ cm-^ for SBS 0335-052E fllzotov et al. Il2006bf ). 



The other rows of Table [TT] hst T^ff, the effective temperature of stellar radiation, and 
Ne, the electron number density, assumed to be constant with radius for all obje cts, except for 
Mrk 9 96, for which we have adopted the number density radial distribution of iThuan et al. 
(120081 ) ■ The fourth row gives /, the filling factor. The remaining input parameters are 
the number ratios of different species to hydrogen, derived mainly from the spectroscopic 
observations in this paper, except for the carbon and silicon abundances. There are no 
strong emission lines of these last two species in the optical and NIR ranges. We have 
therefore calculated their abundances, using the mean re lations of C/Q and Si/Q vs. oxygen 



abundance derived for low-metallicity H ll regions by iGarnett et al. I (1995ajjb|). For the 



ionizing stellar radiation, we have adopted a iKurucz I (jl979[ ) stellar atmosphere model with 



Te//-=50000 K. We have also experimented w ith other forms of the stellar ionizing radiation 



such as those given by Starburst99 models ( jLeitherer et al. 1119991 ) or by models based on 
the Geneva stellar ev olutionary tracks with heavy element mass fractions Z = 0.001 - 0.004 
( iMeynet et al. Ill994j ). In all cases, we were able to reproduce well most of the observed line 
intensities by varying only slightly the input parameters. Thus our results are robust and 
do not depend sensitively on a particular adopted grid of models. 

Most of the input parameters in Table [TTl e.g. Q(H) and chemical composition, are 
constrained by observations and are thus kept unchanged. This is not the case for the filling 
factor / which is therefore varied until the best agreement is found between the observed 
and predicted line intensities. The CLOUDY-predicted fluxes of the emission lines are shown 
in the fifth column of Tables |3] - [71 Comparison of the extinction-corrected observed and 
predicted emission- line fluxes in both the optical and NIR ranges shows that, in general, the 
agreement between the observations and the CLOUDY predictions is satisfactory. This is 
especially true for the strong hydrogen, neutral helium and forbidden lines of doubly ionized 
species. The agreement however is not as good for emission lines of singly ionized and neutral 
species which arise in the outer parts of H ll regions or even beyond the ionized hydrogen 
regions (e.g. the [S ll] emission lines). For example, the line intensities of the [S ll]AA 
6717,6731 emission lines are underpredicted by a factor of 2 for the nearer BCDs (Mrk 71 
and II Zw 40) and by a factor of 3-4 for the more distant ones (SBS 0335-052 and Mrk 930). 
The larger disagreement factor for the more distant BCDs is probably due to the fact that 
their extraction apertures are larger in linear size and include more of the weakly ionized 
and neutral ISM than the apertures of the nearer BCDs. 

Some part of the disagreement can be caused by our too simplistic modeling with the 
CLOUDY code. For all galaxies but Mrk 996, we have adopted a spherically symmetric H 
II region model with a single compact central ionizing source and an uniform density distri- 
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bution. In reality, our galaxies are likely to contain several non- symmetric and non-uniform 
H II regions with different ionization parameters which contribute in different proportions 
to the excitation of high- and low-ionization species in the integrated spectra. With the 
exception of Mrk 996, we do not have enough observational data to produce more complex 
and refined models of the H ll regions. For the sole case o f Mrk 996, we have adop ted a 
non- uniform radial density distribution (the one dervived by iThuan et al. Ill996l . l2008l ). and 
there is indeed a much better agreement between the observed and predicted line intensities, 
particularly for the [S ll]AA 6717,6731 emission lines. 

More sophisticated modeling will certainly improve the overall agreement between ob- 
servations and models. However, these improvements will not change the main conclusion 
that H II region models including only stellar photoionization as an ionizing source can ac- 
count for most of the observed line fluxes, both in the optical and NIR ranges. For the great 
majority of the lines, no additional excitation mechanism such as shocks from stellar winds 
and supernova remnants i s needed. This situ ation is again simi l ar to that in the BCDs II 
Zw 40 and Mrk 59, where IVanzi et al. I mm and llzotov et al. I fco09h found that the ISM 
is mainly photoexcited by stars. However, we note that while some He I emission lines in 
the spectrum of Mrk 996 are reproduced reasonably well, the modeled fluxes of the He I 
A0.707/im and 1.083/im emission lines are several times larger than those observed. These 
two lines are very sensiti ve to coUisiona l and fluorescent excitation. In the very high-density 
H II region of Mrk 996 dThuan et aLlboOsj ). the He I A1.083/im is optically thick. This 
may be a reason why this line has an asymmetric profile (see inset in Fig. H]). Therefore, a 
detailed treatment of the radiative transfer in this line and of the coUisional excitation of He 
I states is needed. 



3.4. The contribution of shock photoionization 

There are other exceptions. One such exception is the He ll AO. 469 /im nebular emission 
line. It is detected in the optical spectra of three BCDs, Mrk 71, Mrk 930 and SBS 0335-052E, 
but cannot be reproduced by the photoionized H ll region models with only stellar radiation 
as the ionizing source. In the two remaining BCDs, II Zw 40 and Mrk 996, this emission 
line is also detected, but it is broad, implying that it is produced by expanding envelopes 
of Wolf-Rayet stars, and not by nebular emission. To reproduce the observed nebular He ll 
AO. 469 nm emission-line flux, a harder ionizing radiation is required. Therefore, for all BCDs 
(except for II Zw 40 and Mrk 996), we have also calculated models of H ll regions with a 
constant number density, which include ionization radiation by fast shocks, in addition to 
the stellar ionizing radiation. We vary the number of ionizing photons produced by fast 
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shocks until the model prediction best matches the observed flux of the nebular He ll AO. 469 
/im emission line. We adopt the spect ral distribution of the ionizing radiation to be that of 
a shock with a velocity of 500 km s~^ ( lAUen et al. II2008I ). For best agreement, we found the 
number of ionizing photons from fast shocks to be respectively 2.5%, 4% and 13% the number 
of ionizing stellar photons for Mrk 71, Mrk 930 and SBS 0335-052E. In the last column of 
Tables m [5] and m we show the emission-line fluxes predicted by combined models including 
both stellar and shock ionizing radiation. These models reproduce well the observed nebular 
He II AO. 469 /im emission line fluxes. They also increase the relative fluxes of the low- 
ionization species, such as [O l] AO. 630 and [S ll] AO. 672, 0.673 fim emission lines, making 
them agree better with the observations. On the other hand, the relative line fluxes of the 
high-ionization species such as [Ne ill] AO. 387, and [S ill] AO. 9533 /im are little changed. This 
is probably because the shock-produced harder ionizing radiation penetrates more efficiently 
into the outer layers of the H ll region, as compared to the softer stellar ionizing radiation. 
This results in more efficient heating of the gas and excitation of the low-ionization species. 



3.5. [Fe ll] line emission 



The [Fe ll] A 1.25 7 and A 1.643 fim emission lines have often been used as indicators of 



the iraportance of shock excit ation relative to that of photoionization (IMoorwood &: Oliva 



1988; 


Oliva et al. 


1990, 


2001) 



1990L |2001[ ). All of our galaxies either show one of them or both. In the 
Galaxy, the intensities of these lines are considerably smaller in regions where photoionization 
dominates than in those where shock excitation dominates, such as in SN remnants. For 
example, the ratio of the [Fe ll] A1.643 /ii n line to Br7 is 0.06 in the Orion nebula, whi le it 



is more than 30 in Galactic SN remnants (IMoorwood fc Oliva Ill988l : lOliva et al. Ill990l ). 



The [Fe ll] A1.257 fim and A1.643 /im lines to Br7 ratios have the very low values of 
0.05 - 0.1 in II Zw 40, Mrk 71 and SBS 0335-052E, similar to the value of 0.055 found 



by llzotov et al. I (120091 ) in Mrk 59, and to the one in the Orion nebula. However, in the 
two other galaxies, Mrk 930 and Mrk 996, which have lower excitation H ll regions, the [Fe 
Il]/Br7 ratios are considerably higher, with values ~ 0.5, but still much lower than the value 
of >30 in Galactic SN remnants. Does this mean that SN shock excitation does not play a 
role in all our BCDs? 



Examination of Table fTH shows that, even without any contribution of shock excitation 
from SNe, the stellar photoionization CLOUDY models overpredict the [Fe ll] A1.257 and 
A1.643 /im line fluxes by factors of ~1.5 - 3. An exception is Mrk 930, where the extinction- 
corrected [Fe ll] line fluxes are a factor of ~ 2 greater than the modeled ones. If shock 
ionizing radiation is included (last column of Table [T^ . then the modeled [Fe ll] emission 
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line fluxes are more comparable to the observed ones in Mrk 930. However, they are much 
higher than the observed fluxes in other galaxies. It appears that, for H ll regions in II Zw 
40, Mrk 71, Mrk 996 and SBS 0335-052E, there is no need to invoke any shock ionizing 
radiation. However, shocks may play an important role in the excitation of [Fe ll] emission 
lines in Mrk 930. 

The overprediction of the [Fe ll] fluxes in the four BCDs other than Mrk 930 may be due 
to aperture effects. The small apertures used in the optical and NIR observations (the sht 
widths are respectively 1'.'5 and I'.'l) include only the inner parts of the H ll region while singly 
ionized iron is mainly located in its outer parts. This conclusion appears to be supported 
by the fact that the discrepancy factor is largest for the nearest BCD, Mrk 71, at a distance 
of 3.4 Mpc, while it is lower for more distant galaxies, with distances larger than 10 Mpc . 



The aperture effect is the smallest for the most distant BCD, Mrk 930. Ilzotov et al. I (|2009[ ) 
have proposed another explanation for the BCD Mrk 59, at a distance of 10.7 Mpc. They 
argued that the overprediction by a factor of ~ 2 of the [Fe ll] 1.257 and 1.643 /im emission 



( Oliva et al. 


2001; 


Izotov et al. 


2006a) 



depletion process is likely to be most efficient in the outer zones of the H ll region where 
the [Fe ll] emission lines arise. These have less intense ionizing radiation as compared to 
the inner zones where the ionizing radiation is more intense and where the [Fe ill] optical 
emission lines arise. It is likely that both aperture effects and depletion of iron onto dust 
play a role in accounting for the discrepancies between the modeled and observed [Fe ll] 
fluxes. 



3.6. Comparison with mid-infrared observations: MIR extinction 



We have shown above that the observed relative line flux ratios in both the optical and 
NIR ranges in all our objects can be satisfactorily accounted for with the same extinction 
coefficient C(H/3) [or the same extinction y4(y)]. In other words, there is no more hidden 
star formation in the NIR range as compared to in the optical range in regions traced by 
emission lines, i.e. in ionized gas regions. A question then arises: would we see more hidden 
star formation at longer wavelengths, in the mid-infrared (MIR) range? 

To investigate the extir iction in th e MIR ran ge, we have used the pu blished Spitzer 



M IR emission line fluxes of IWu et al. I ( 120061 ) and lO'Halloran et a. 



of O'Halloran et a 



((20061) for Mrk 930, of iThuan et al. I (120081 ) for Mrk 996, and of 



(120061 ) for II Zw 40, 



Houck et al. I (|2004[ ) for SBS 0335-052E. Spitzer observations of Mrk 71 exist but have not yet 



been published. Since the Spitzer spectra were obtained within significantly larger apertures 
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as compared to the optical and NIR spectra, we have normahzed the MIR emission-hne 
fluxes to the optical and NIR fluxes by multiplying them by the ratios of the extinction- 
corrected optical H/5 fluxes listed in Tables 3-7 to those obtained from large aperture spec- 
troscopic observations or from ph otometric observati ons. Thus, for II Zw 40 and Mrk 93 
we have adopted the H/3 fluxes of iLagos et al. I (120071 ) and iMoustakas fc Kennicutt I (120061 ) . 
respectively. These fluxes were then corrected for extinction with the extinction coefficient 
C(II/3) derived in this paper from the optical spectra. For Mrk 996 and SBS 335-0 52E, 



we have adop t ed resp ectively the extinction-corrected H/3 fluxes of iThuan et al. I (120081 ) and 
Izotov et al. I (l2006al ). The multiplicative factors are respectively 0.32, 0.27, 0.31, and 0.50 



for II Zw 40, Mrk 930, Mrk 996 and SBS 0335-052E. No extinction correction has been 
applied to the MIR data as it is presumably small in this wavelength range. 



The normalized MIR emission line fluxes are shown in Table [13] together with the flux 
predictions by the CLOUDY models discussed in Section [373] for the optical and NIR ranges. 
Comparison of the MIR observed fluxes with those predicted by the pure stellar ionizing 
radiation model shows agreement to within a factor of 2. The disagreement is the worst for 
the low ionization lines [Ne ll] A12.81 fxm in II Zw 40 and [Fe ll] A25.99 fxm in II Zw 40 
and Mrk 930, the fluxes of which are predicted to be considerably smaller than the observed 
ones. Furthermore, the presence of the high-ionization emission line [O iv] A25.89 /xm in 



Mrk 996 cannot be accounted for by a pure stellar ionizing radiation model (IThuan et al 



2008]). The high observed flux of [Fe ll] A25.99 /im in Mrk 930 cannot be explained by shock 
excitation, but it can be understood if we note that neutral iron has an ionization potential 
smaller than that of hydrogen. Therefore, singly ionized iron may exist in weakly ionized 
cool interstellar medium, outside the H ll regions, and contribute to the [Fe ll] A25.99 fim 
emission, a contribution which is not accounted for in CLOUDY mod e ls. A s for the high- 
ionization [O iv] A25.89 fxm emission line in Mrk 996, Thuan et al. ( 20081 ) found that a 
CLOUDY model where the fraction of ionizing photons from shocks is only 3% fits well all 
the MIR emission- lines of Mrk 996, including [O iv] A25.89 /zm. We give the predictions 
of that model for Mrk 996 in Table [131 The agreement within a factor of 2 between the 
observed fluxes of the MIR emission-lines and the predicted ones of CLOUDY models, based 
on the extinction-corrected optical emission lines, implies that there is not considerably more 
hidden star formation seen in the MIR range, as compared to in the optical and NIR ranges. 
This implies that, in low-metallicity BCDs, the H ll regions with emission lines detected in 
the MIR range are in general relatively transparent, so they are also seen in the optical and 
NIR ranges. 

At first glan c e, thi s conc lusion appears to be in flagrant disagreement with the findings 



of [Thuan et al. I (11999] ) and ]Houck et al. I (12004] ) who inferred highly obscured regions in 
SBS 0335-052E, with an A{V) ~ 15 mag, from modeling of the MIR continuum emission 
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produced by warm dust. Based on these observations, iThuan et al. I ( 119991 ) concluded that 
some 75% of the star formation in SBS 0335-052E is hidden. These two apparently discrepant 
results can be reconciled if the region where resides the warm dust responsible for the MIR 
emission is different from the ionized ga s region where the emissi on lines arise. Studies of 
the H II regions in the Milky Way (e.g., iLebouteiller et al. 1120071 ) show a good correlation 
between e.g. the [S iv] A10.51/xm emission line flux and the warm dust continuum, implying 
that the warm dust resides in the ionized medium. However, the physical conditions in SBS 
0335-052E are very different from those in the Milky Way H ll regions. In particular, star 
formation in SBS 0335-052E occurs in several cq r apact SSCs containiri g several thousand 



low-metallicity O stars (e.g., iThuan et al. 1 119971 : I Johnson et al. II2009I ). responsible for a 



very high UV flux in the H II regions around the SSCs. It is doubtful whether dust grains 
can survive in such an extreme environment. Furthermore, if the warm dust is in the same 
region as the ionized gas, we would expect an excess of the MIR emission line fluxes as 
compared to the predicted ones. No such excess is found, at least for the [Ne ill] A15.55 
/im emission line (Table [T3|l . The dust is likely located in the neutral gas that surrounds 
the compact H ll regions and which does not contribute to the observed MIR emission lines. 
In other words, the MIR emission lines and continuum arise in different regions, the first 
ones in relatively transparent H ll regions and the second one in highly extincted regions 
surrounding the H II regions. 



4. SUMMARY AND CONCLUSIONS 

We present here near-infrared (NIR) spectroscopic observations of five blue compact 
dwarf (BCD) galaxies, II Zw 40, Mrk 71, Mrk 930, Mrk 996 and SBS 0335-052E, obtained 
with the 3.5m APO telescope. The NIR data have been supplemented by optical spectra 
obtained with the 3.5m APO and the 2.1m KPNO telescopes, resulting in a total wavelength 
coverage of 0.36 - 2.46 /xm. The overlap in the 0.90 - 0.97 /im wavelength range allows to 
adjust the continuum levels of the NIR spectra to the optical ones by insuring that the flux 
of the [S III] AO. 953 /im emission line is the same in both types of spectra. 

We have arrived at the following conclusions: 

1. In all BCDs, the extinction coefficient C(H/3) = 0A6xA{V) derived in the optical is 
the same as the one derived in the NIR for the ionized gas regions traced by emission lines. 
The NIR emission lines do not probe more extinct regions with hidden star formation as 
compared to the optical emission lines. Using Spitzer MIR data, we have also found that the 
extinction A{V) inferred from MIR observations is similar to the one derived from optical 
and NIR observations. Thus, the emission-line spectra of low-metallicity BCDs in the whole 



-17- 



~ 0.36 - 25 iim wavelength range originates in relatively transparent H ll regions. The large 
extinction deduced from the MIR continuum emission in the BCD SBS 0335-052E probably 
arises from neutral regions surrounding the H ll regions. 

2. We have detected molecular hydrogen emission hues in the NIR spectra of all BCDs. 
Comparison of the observed fluxes with modeled ones suggests that the main excitation 
mechanism of H2 emission in all BCDs is fluorescence. 

3. A CLOUDY model with a pure stellar ionizing radiation reproduces well the observed 
fluxes of most emission hues in both optical and NIR ranges. Shocks are hkely present in 
the H II region, but play a minor role in the ionization, at the level of 10% or less, except 
for high-ionization emission lines, such as He ll A0.469 /im in the optical range and [O iv] 
A25.89 /xm in the MIR range. 

4. CLOUDY stellar and shock ionizing models show that the [Fe ll] A1.257 and A1.643 

/xm emission lines, often used as SN shock indicators in low-excitation high-metallicity star- 
burst galaxies, cannot play such a role in high-excitation low-metallicity H ll regions, such 
as those in the BCDs II Zw 40, Mrk 71, Mrk 996 and SBS 0335-052E. However, we find that 
shocks could be important contributors to the emission of [Fe ll] lines in Mrk 930. 

The 3.5 APO time was available thanks to a grant from the Frank Levinson Fund of 
the Silicon Valley Community Foundation to the Astronomy Department of the University 
of Virginia. Y.I.I, is grateful to the staff of the Astronomy Department at the University of 
Virginia for their warm hospitality. T.X.T acknowledges the support of NASA. 
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Table 1. General characteristics of galaxies 



Object 


R.A. (J2000.0) 


Dec. (J2000.0) 


Redshift 


Distance*^ 


12+logO/H^ 


II Zw 40 


05 55 42.6 


+03 23 32 


0.00263 


10.4 


8.00 


Mrk 71 


07 28 42.8 


+69 11 21 


0.00036 


3.4 


7.85 


Mrk 930 


23 31 58.3 


+28 56 50 


0.01830 


77.4 


8.08 


Mrk 996 


01 27 35.5 


-06 19 36 


0.00541 


21.8 


8.10 


SBS 0335-052E 


03 37 44.0 


-05 02 40 


0.01348 


53.8 


7.29 



^Distances of all galaxies but Mrk 71 are obtained from the radial velocities corrected 
for the Virgo infall and are tak en from the NASA/IPA C Extragalactic Database (NED). 
The distance of Mrk 71 is from iThuan &: Izotov I (|2005l ). 

'^this paper. 



Table 2. Journal of observations 



Object 




Date 


Exposure (s) 


Extraction 
aperture^ 


Airmass 








Near-infrared observations 








II Zw 40 


25 


Nov. 


2009 


2000 (4x300 + 4x200) 


l'.'lx6" 


1 


7 


Mrk 71 


25 


Nov. 


2009 


1200 (4x300) 


l'.'lx6" 


1 


8 


Mrk 930 


25 


Nov. 


2009 


1200 (4x300) 


l'.'lx6" 


1 





Mrk 996 


25 


Nov. 


2009 


2400 (8x300) 


l'.'lx6" 


1 


7 


SBS 0335-052E 


25 


Nov. 


2009 


3600 (12x300) 


l'.'lx6" 


1 


7 








Optical observations 








II Zw 40 


07 


Feb. 


2010 


1200 


l'.'5x4" 


1 


2 


Mrk 71 


20 


Nov. 


2009 


660 


l'.'5x4" 


1 


3 


Mrk 930 


16 


Dec. 


1996 


3600 


2"x6'.'9 


1 





Mrk 996 


06 


Nov. 


2008 


1200 


l'.'5x4" 


1 


7 


SBS 0335-052E 


20 


Nov. 


2009 


1200 


l'.'5x4" 


1 


3 



Extraction aperture for the one-dimensional spectrum. 



Table 3. Fluxes of optical and NIR emission lines in II Zw 40 



CLOUDY 



Ion 


100xF(A)/F(H^)^ 


100x/(A)//(H/3)'^ 


100x/(A)//(H^)^ 


stars 




a) Optical spectrum 






0.373 [0 ii] 


23.87ib0.66 


59.06=b2.03 


55.59zb2.03 


53.28 


0.380 HIO 


3.30±0.39 


7.53ibl.29 


7.49zbl.29 


5.68 


0.384 H9 


3.91±0.34 


8.70ibl.27 


8.67zbl.28 


7.70 


0.387 [Ne iii] 


28.90±0.56 


61.52ibl.94 


62.24zbl.94 


60.73 


0.389 He i+HS 


8.52±0.34 


18.14ibl.39 


18.22zbl.39 


17.84 


0.389 [Ne iii]+H7 


17.70ib0.41 


35.54zbl.56 


35.89zbl.57 


34.62 


0.410 RS 


17.89±0.37 


or* ■! 1 -t 0/^ 

32.61zbl.36 


32.04zbl.33 


26.34 


0.434 H7 


32.89±0.54 


A r\ A / \ 1-1 A 

49.49zbl.48 


48.92zbl.47 


47.42 


0.436 [0 III] 


8.67±0.22 


12.77zb0.70 


13.04±0.71 


11.34 


0.447 He I 


2.83±0.17 


3.81zb0.38 


3.90zb0.38 


3.95 


0.471 [Ar iv]+He I 


2.68±0.14 


a r\ci 1 n or\ 

2.98±0.30 


3.07±0.30 


2.77 


0.474 [Ar iv] 


1.18±0.13 


1.29±0.20 


1.33zb0.21 


1.69 


0.486 H/3 


1 r\f\ c\c\ 11 AC\ 

100.00±1.49 


100.00±1.59 


100.00zbl.60 


100.00 


0.496 [0 III] 


258.52±3.79 


241.50zb3.72 


249.58zb3.84 


226.58 


0.501 [0 III] 


797.24±11.6 


721.72zbl0.6 


746.46±11.0 


682.00 


0.588 He I 


13.75ib0.24 


8.06zb0.35 


8.43zb0.36 


9.99 


0.630 [0 i] 


z.4U±U.10 


1.20±0.13 


1.26zb0.13 


0.68 


0.631 [S III] 


2.98zb0.11 


1.49zb0.14 


1.56zb0.14 


1.70 


0.636 [0 i] 


l.MzbO.lO 


0.56±0.09 


0.59zb0.09 


0.22 


0.656 Ha 


622.51=b9.07 


283.94zb4.56 


283.94zb4.56 


282.41 


0.658 [N 11] 


9.54ib0.17 


4.32zb0.23 


4.55zb0.23 


2.33 


0.668 He I 


7.26ib0.14 


3.18zb0.19 


3.35zb0.20 


2.72 


0.672 [S 11] 


10.72ib0.19 


4.63±0.23 


4.88zb0.24 


2.64 


0.673 [S 11] 


10.20ib0.18 


4.39zb0.22 


4.62zb0.23 


2.00 


0.707 He I 


13.13ib0.21 


5.01zb0.23 


5.30zb0.24 


4.20 


0.714 [Ar III] 


21.55ib0.34 


8.03zb0.29 


8.49zb0.30 


8.91 


0.728 He i 


2.18ib0.09 


0.77zb0.08 


0.82zb0.09 


0.68 


0.732 [0 11] 


4.33ib0.11 


1.51zb0.12 


1.60zb0.12 


0.98 


0.733 [0 11] 


3.79±0.14 


1.32zb0.11 


1.40zb0.12 


0.78 


0.775 [Ar III] 


6.20±0.15 


1.86zb0.12 


1.98±0.13 


2.15 


0.805 [CI IV] 


1.26±0.08 


0.34±0.05 


0.37±0.05 


0.39 
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Table 3 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H^)^ 


100x/(A)/I(H/3)'^ 


100xI(A)//(H^)^ 


stars 


0.847 Pal7 


1.32±0.09 


0.33±0.09 


0.35±0.09 


0.29 


0.850 Palo 


2.30±0.11 


0.56±0.09 


0.59±0.10 


0.35 


0.855 Pal5 


2.19±0.10 


0.53±0.09 


0.56±0.09 


0.45 


0.860 Pal4 


3.25±0.11 


0.77±0.09 


0.81±0.10 


0.98 


0.867 Pais 


3.02±0.11 


0.70±0.09 


0.74±0.10 


1.02 


0.875 Pal2 


4.76±0.13 


1.07±0.10 


1.14±0.10 


1.20 


0.886 Pall 


5.11±0.15 


1.12±0.10 


1.19=b0.10 


1.48 


0.901 PalO 


9.72±0.21 


2.05±0.12 


2.19±0.12 


1.92 


0.907 [S III] 


62.59±0.95 


12.97±0.33 


13.93±0.35 


13.30 


0.923 Pa9 


10.48±0.23 


2.12±0.12 


2.27±0.13 


2.58 


0.953 [S III] 


141.62±2.14 


26.92±0.59 


28.97±0.63 


32.98 


0.955 Pae 


19.30±0.50 


3.67±0.16 


3.94±0.17 


3.65 




b) Near infrared spectrum 






0.907 [S III] 


59.60±8.72 


12.34±1.81 


13.26±1.95 


13.30 


0.953 [S III] 


141.62±7.73 


26.92±1.52 


28.97±1.63 


32.98 


0.955 Pae 


16.26±0.96 


3.09±0.20 


3.32±0.22 


3.65 


1.003 He I 


1.84±0.17 


0.32±0.04 


0.35±0.04 


0.12 


1.005 Pa(5 


26.89±1.47 


4.71±0.28 


5.07=b0.30 


5.45 


1.029 [S ii] 


0.55±0.10 


0.09±0.02 


0.10±0.02 


0.08 


1.032 [S ll] 


0.49±0.06 


0.08±0.02 


0.09±0.02 


0.11 


1.083 He I 


190.47±10.4 


29.93±1.69 


32.39±1.83 


26.94 


1.091 He I 


2.75±0.17 


0.43±0.04 


0.46±0.05 


0.28 


1.094 Pa7 


47.21±2.57 


7.33±0.42 


7.92±0.45 


8.74 


1.197 He I 


1.74±0.11 


0.24±0.03 


0.26±0.03 


0.19 


1.207 H2 3-1 S(2) 


O.MibO.Ol 


0.02±0.01 


0.02±0.01 




1.233 H2 3-1 S(l) 


0.30±0.04 


0.04±0.02 


0.04±0.02 




1.238 H2 2-0 Q(l) 


0.43±0.06 


0.06±0.03 


0.06±0.03 




1.253 He I 


1.57±0.10 


0.21±0.03 


0.23±0.03 


0.25 


1.257 [Fe 11] 


2.13±0.13 


0.28±0.03 


0.30±0.03 


0.39 


1.279 He I 


6.68±0.37 


0.86±0.06 


0.94±0.07 


0.69 


1.282 Pa/3 


94.07±5.12 


12.13±0.69 


13.19±0.75 


15.55 


1.311 H2 4-2 S(l) 


0.28±0.04 


0.04±0.02 


0.04±0.02 
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Table 3 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H^)^ 


100x/(A)/I(H/3)'^ 


100xI(A)//(H^)^ 


stars 


1.314 H2 3-1 Q(l) 


0.31±0.01 


0.04±0.02 


0.04±0.02 




1.500 Br24 


0.72±0.06 


0.08±0.02 


0.09±0.02 


0.20 


1.504 Br23 


0.79±0.06 


0.09±0.02 


0.10±0.02 


0.22 


1.509 Br22 


1.12ib0.07 


0.13±0.02 


0.14±0.03 


0.23 


1.514 Br21 


1.18±0.07 


0.14±0.03 


0.15±0.03 


0.09 


1.519 Br20 


1.45±0.08 


0.17±0.03 


0.18ib0.03 


0.10 


1.526 Brl9 


0.98±0.06 


0.11±0.02 


0.12±0.02 


0.11 


1.535 Brl8 


1.56±0.09 


0.18±0.03 


yA -1 yA 1 

0.19±0.03 


0.12 


1.544 Brl7 


1.53±0.09 


0.17±0.03 


0.18±0.03 


0.14 


1.556 Brl6 


2.25=b0.13 


0.25±0.03 


0.27±0.03 


0.17 


1.570 Brl5 


2.18ib0.12 


0.24±0.03 


0.26±0.03 


0.22 


1.588 Brl4 


2.63±0.15 


0.29±0.03 


0.31±0.03 


0.47 


1.601 H2 6-4 Q(l) 


0.24±0.03 


0.03±0.01 


yA y^ 1 -1 

0.03±0.01 




1.611 Brl3 


3.73±0.21 


0.40±0.04 


0.43±0.04 


0.49 


1.641 Brl2 


4.25±0.23 


0.45±0.04 


0.49±0.04 


0.57 


1.643 [Fe iij 


2.05±0.11 


0.22±0.03 


0.24±0.03 


0.38 


1.681 Brll 


6.10±0.33 


0.63±0.05 


0.69±0.05 


0.71 


1.700 He I 


2.77±0.15 


0.28=b0.03 


0.31±0.03 


0.28 


1.736 Brio 


8.21±0.45 


0.83=b0.06 


y^ yA 1 yA jA 

0.91±0.06 


0.90 


1.818 Br9 


14.10±0.77 


1.39±0.09 


1.52±0.10 


1.24 


1.944 Bro 


13.26±0.72 


1.26±0.08 


1.38±0.09 


1.75 


1.957 H2 1-0 S(3) 


0.69±0.05 


0.07±0.02 


0.07±0.02 




2.034 H2 1-0 S(2) 


0.38±0.03 


0.04±0.01 


0.04±0.01 




2.058 He i 


7.27±0.40 


0.67±0.05 


0.73±0.05 


1.45 


2.073 H2 2-1 S(3) 


0.33±0.02 


0.03±0.01 


0.03±0.01 




2.113 He I 


0.72±0.04 


0.07±0.01 


0.07±0.01 


0.08 


2.122 H2 1-0 S(l) 


1.13±0.06 


0.10±0.02 


0.11±0.02 




2.165 Br7 


27.90±1.52 


2.51±0.15 


2.75±0.16 


2.61 


2.218 [Fe III] 


0.55±0.03 


0.05±0.01 


0.05±0.01 




2.223 H2 1-0 S(0) 


0.40±0.02 


0.04±0.01 


0.05±0.01 




2.248 H2 2-1 S(l) 


0.61±0.05 


0.06±0.02 


0.06±0.02 




2.355 H2 2-1 S(0) 


0.22±0.02 


0.02±0.01 


0.02±0.01 





Table 3 — Continued 



CLOUDY 

Ion 100xF(A)/F(H^)^ 100x/(A)//(H;S)'^ 100x/(A)//(H^)^ sta^s 



^F(H/3)=(1.097±0.004)xlO-i^ erg s"^ cm-^. 

''Corrected for extinction with C(H/5)=1.165 and EW(abs)=0.250A. 
'^Corrected for extinction with C(H;S)=1.135 and EW(abs)=0.100A. 
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Table 4. Fluxes of optical and NIR emission lines in Mrk 71 



CLOUDY 



Ion 


W0xF{X)/F{Rl3)'' 


100x/(A)//(H^)'^ 


100x/(A)//(H/3)^ 


stars 


stars+shock 




a) Optical spectrum 








0.373 [0 ii] 


37.61±0.56 


A A nr* 11 11 

44.09±1.11 


44.53zfcl.ll 


35.52 


36.54 


0.375 H12 


2.58±0.11 


3.35ifc0.30 


3.27zfc0.29 


3.54 


3.53 


0.377 Hll 


3.02±0.11 


3.84ifc0.32 


3.77zfc0.31 


4.38 


4.37 


0.380 HIO 


3.98±0.11 


4.95ifc0.35 


4.89zfc0.35 


5.68 


5.67 


0.382 He i 


0.79±0.09 


0.91ifc0.16 


0.92zfc0.16 


1.04 


1.03 


0.384 H9 


6.06±0.13 


7.35=fc0.43 


7.32zfc0.42 


7.69 


7.67 


0.387 [Ne iii] 


52.45rfc0.76 


60.55±1.32 


61.27zfcl.32 


60.40 


60.47 


0.389 He I+H8 


17.94ifc0.28 


21.00±0.73 


21.15zfc0.72 


17.52 


17.49 


0.389 [Ne iil]+H7 


27.94ifc0.41 


32.16±0.91 


32.49zfc0.90 


34.51 


34.53 


0.403 He I 


1.67ifc0.08 


1.89±0.21 


1.92zfc0.21 


1.90 


1.89 


0.407 [S ii] 


0.59±0.10 


0.67±0.14 


0.68zfc0.14 


0.26 


0.32 


0.410 B.S 


24.13±0.35 


27.28±0.82 


26.80zfc0.80 


26.34 


26.29 


0.434 H7 


A f Ar\ \ r\ /^f 

45.40±0.65 


Af\ cxr \ -A 1 -| 

49.25±1.11 


48.69±1.10 


47.45 


47.42 


0.436 [0 III] 


15.18ifc0.22 


16.32zfc0.61 


16.70zfc0.61 


12.63 


12.88 


0.439 He i 


0.49ifc0.04 


0.52ifc0.11 


0.53±0.11 


0.51 


0.51 


0.447 He i 


3.76ifc0.07 


3.98ifc0.29 


4.07zfc0.30 


3.95 


3.92 


0.469 He li 


0.46±0.04 


0.47±0.10 


0.49zfc0.10 


0.00 


0.53 


0.471 [Ar lv]+He I 


3.22zfc0.06 


3.28zfc0.26 


3.38zfc0.27 


2.98 


3.01 


0.474 [Ar iv] 


2.26zfc0.05 


2.29=fc0.22 


2.36zfc0.22 


1.83 


1.87 


0.486 H/3 


100.00ifcl.42 


100.00zfcl.43 


100.00zfcl.43 


100.00 


100.00 


0.492 He i 


0.77ifc0.03 


0.76ifc0.12 


0.79zfc0.13 


1.05 


1.04 


0.496 [0 III] 


243.93±3.46 


240.18±3.56 


248.25zfc3.66 


229.19 


229.53 


0.501 [0 III] 


742.38ifcl0.5 


726.55±10.4 


751.56±10.7 


689.86 


690.90 


0.575 [N 11] 


0.23±0.01 


0.21±0.06 


0.22zfc0.06 


0.02 


0.02 


0.588 He I 


12.59ifc0.18 


ll.33zfc0.46 


ll.85zfc0.48 


10.02 


9.94 


0.630 [0 i] 


1.05ifc0.02 


0.93zfc0.13 


0.96±0.13 


0.45 


0.68 


0.631 [S III] 


1.55ifc0.02 


1.36zfc0.16 


1.42zfc0.16 


1.44 


1.43 


0.636 [0 l] 


0.37ifc0.02 


0.32zfc0.08 


0.34zfc0.08 


0.15 


0.22 


0.656 Ha 


324.850ifc4.61 


278.68zfc4.41 


278.68zfc4.41 


281.65 


281.90 


0.658 [N 11] 


1.47±0.03 


1.26zfc0.15 


1.33±0.15 


0.62 


0.66 


0.668 He i 


3.61±0.05 


3.07±0.23 


3.24zfc0.24 


2.70 


2.68 
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Table 4 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H^)'^ 


100xI(A)//(H^)'' 


100x/(A)/I(H/3)^ 


stars 


stars+stiock 


0.672 [S ii] 


4.14±0.06 


3.51±0.25 


3.70±0.26 


1.79 


1.96 


0.673 [S ii] 


3.43±0.05 


2.91±0.23 


3.07±0.23 


1.36 


1.48 


0.707 He I 


4.59±0.07 


3.80±0.26 


4.02±0.27 


4.80 


4.73 


0.714 [Ar illj 


6.54±0.10 


5.39±0.31 


5.71±0.32 


6.49 


6.48 


0.728 He i 


0.74±0.02 


0.61±0.10 


0.64±0.11 


0.69 


0.68 


0.732 [0 ii] 


1.05±0.02 


0.85ib0.12 


0.90±0.12 


0.68 


0.70 


0.733 [0 ii] 


0.93±0.01 


0.76ib0.11 


0.80±0.12 


0.54 


0.56 


0.775 [Ar iii] 


1.62±0.03 


1.28±0.15 


1.36±0.15 


1.57 


1.57 


0.805 [CI iv] 


0.48±0.02 


0.37±0.08 


0.40±0.08 


0.38 


0.38 


0.835 Pa23 


0.21±0.02 


0.20±0.06 


0.20±0.06 


0.45 


0.45 


0.836 Pa22 


0.29±0.02 


0.26±0.06 


0.26±0.07 


0.48 


0.48 


0.837 Pa21 


0.24±0.02 


0.22±0.06 


0.22±0.06 


0.19 


0.19 


0.839 Pa20 


0.27±0.02 


0.23ib0.06 


0.24ib0.06 


0.20 


0.20 


0.841 Pal9 


0.34±0.02 


0.29±0.07 


0.30±0.07 


0.22 


0.22 


0.844 Pais 


0.39±0.02 


0.33±0.07 


0.35±0.07 


0.25 


0.25 


0.845 I 


0.22±0.01 


0.17±0.05 


0.18±0.05 


0.00 


0.00 


0.847 Pal7 


0.42±0.02 


0.35±0.07 


0.37±0.08 


0.29 


0.29 


0.850 Pal6 


0.66±0.02 


0.53±0.09 


0.56ib0.10 


0.35 


0.35 


0.855 Pal5 


0.68±0.02 


0.54±0.09 


0.57±0.10 


0.45 


0.45 


0.860 Pal4 


0.83±0.02 


0.65±0.10 


0.69±0.11 


0.98 


0.98 


0.867 Pal3 


0.99±0.02 


0.78±0.11 


0.82±0.12 


1.02 


1.02 


0.875 Pal2 


1.33±0.03 


1.02±0.13 


1.08±0.13 


1.19 


1.19 


0.886 Pall 


1.78±0.03 


1.35±0.14 


1.43±0.15 


1.48 


1.48 


0.901 PalO 


1.57±0.03 


1.19±0.14 


1.27±0.14 


1.91 


1.91 


0.907 [S III] 


11.63±0.17 


8.59±0.38 


9.22±0.40 


10.46 


10.47 


0.923 Pa9 


3.17±0.06 


2.34±0.19 


2.51±0.20 


2.57 


2.56 


0.953 [S III] 


31.61±0.45 


22.91±0.67 


24.64±0.71 


25.94 


25.95 


0.955 Pae 


5.06±0.07 


3.68±0.24 


3.96±0.25 


3.63 


3.62 




b) Near infrared spectrum 








0.907 [S III] 


15.88±1.46 


11.69±1.09 


12.56±1.17 


10.46 


10.47 


0.953 [S III] 


31.61±1.72 


22.91±1.30 


24.64±1.40 


25.94 


25.95 


0.955 Pae 


4.68±0.27 


3.42±0.26 


3.67±0.27 


3.63 


3.62 



Table 4 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H^)'^ 


100x/(A)//(H/3)'^ 


100xI(A)//(H^)'= 


stars 


stars-bshock 


1.003 He I 


0.23=b0.02 


0.17ib0.05 


0.18ib0.05 


0.12 


0.12 


1.005 Pa6 


7.05=b0.39 


/A 4 1 /A lA fA 

5.04ib0.33 


5.44ib0.36 


5.42 


5.41 


1.032 [S ii] 


0.11±0.02 


0.08ib0.03 


0.08ib0.04 


0.08 


0.09 


1.083 He I 


45.34±2.46 


31.64ibl.79 


A t"! J 1 -1 / \ A 

34.24ibl.94 


27.71 


27.53 


1.091 He I 


0.51±0.03 


0.36ib0.07 


0.39ib0.08 


0.28 


0.28 


1.094 Pa7 


12.47±0.68 


8.70ib0.52 


9.41ib0.56 


8.69 


8.67 


1.162 H2 2-0 S(l) 


0.03=b0.01 


0.02ib0.01 


0.02ib0.01 






1.197 He I 


0.29±0.02 


/A 1 yA vA 

0.20ib0.05 


0.21ib0.06 


0.19 


0.19 


1.207 H2 3-1 S(2) 


0.06±0.01 


0.04ib0.02 


0.04ib0.02 






1.233 H2 3-1 S(l) 


0.12±0.01 


0.08±0.03 


0.09ib0.03 






1.238 H2 2-0 Q(l) 


O.llibO.Ol 


0.07ib0.03 


0.08ib0.03 






1.253 He I 


0.21±0.02 


0.14ib0.05 


0.16ib0.05 


0.28 


0.28 


1.257 [Fe 11] 


VA £^ 1 iA /"X -J 

0.12±0.01 


0.08ib0.03 


vA yA/A 1 /A iA .J 

0.09ib0.04 


0.28 


0.41 


1.279 He I 


1.19±0.07 


O.SOibO.ll 


0.87ib0.11 


0.69 


0.68 


1.282 Pa/3 


21.73±1.18 


14.61ib0.85 


15.88ib0.92 


15.45 


15.43 


1.311 H2 4-2 S(l) 


O.llibO.Ol 


0.07±0.03 


0.08ib0.03 






1.314 H2 3-1 Q(l) 


O.llibO.Ol 


0.07ib0.03 


0.08ib0.03 






1.514 Br21 


0.15ib0.01 


iA -* VA 1 /'\ A 

0.10ib0.04 


0.11ib0.04 


0.09 


0.09 


1.519 Br20 


0.17ib0.01 


0.12ib0.04 


0.13ib0.04 


0.10 


0.10 


1.526 Brl9 


0.25ib0.02 


0.17ib0.05 


0.18ib0.05 


0.10 


0.10 


1.535 Brl8 


0.21ib0.01 


0.14±0.04 


0.16ib0.05 


0.12 


0.12 


1.544 Brl7 


0.24ib0.02 


0.17ib0.05 


0.18ib0.05 


0.14 


0.14 


1.556 Brl6 


0.37ib0.02 


0.25ib0.06 


0.27ib0.06 


0.17 


0.17 


1.570 Brl5 


0.46±0.03 


0.31±0.06 


0.33±0.07 


0.22 


0.22 


1.588 Brl4 


0.44ib0.03 


0.29±0.06 


0.32±0.07 


0.47 


0.47 


1.601 H2 6-4 Q(l) 


0.03ib0.01 


0.02±0.01 


0.02ib0.01 






1.611 Brl3 


0.61ib0.03 


0.40±0.07 


0.43ib0.08 


0.49 


0.49 


1.641 Brl2 


0.69ib0.04 


0.45ib0.08 


0.49ib0.08 


0.57 


0.57 


1.681 Brll 


0.95±0.05 


0.61±0.09 


0.67±0.10 


0.71 


0.71 


1.700 He I 


0.22ib0.01 


0.14ib0.04 


0.16±0.05 


0.28 


0.28 


1.736 Brio 


1.32±0.07 


0.85±0.11 


0.93±0.12 


0.91 


0.91 


1.748 He I 


O.lOiO.Ol 


0.06±0.03 


0.07±0.03 


0.00 


0.00 
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Table 4 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H/3)'^ 


100xI(A)//(H/3)^ 


100xI(A)//(H^)^ 


stars 


stars-|-shock 


1.818 Br9 


2.19±0.12 


1.40±0.14 


1.54±0.15 


1.23 


1.23 


1.944 Br 5 


2.72±0.15 


1.72±0.16 


1.89±0.17 


1.74 


1.73 


1.957 H2 1-0 S(3) 


0.30±0.03 


0.19±0.05 


0.21±0.05 






2.034 H2 1-0 S(2) 


0.07±0.01 


0.05±0.02 


0.05±0.03 






2.058 He i 


1.39±0.08 


0.88±0.11 


0.96±0.12 


1.45 


1.44 


2.073 H2 2-1 S(3) 


0.07±0.01 


0.05±0.02 


0.05±0.03 






2.122 H2 1-0 S(l) 


0.33±0.02 


0.21±0.05 


0.23±0.06 






2.165 Br7 


3.51±0.19 


2.20±0.18 


2.42±0.20 


2.59 


2.58 


2.223 H2 1-0 S(0) 


0.16±0.01 


0.10±0.04 


0.11±0.04 






2.248 H2 2-1 S(l) 


0.15±0.01 


0.10±0.04 


0.11±0.04 






2.355 H2 2-1 S(0) 


0.08±0.01 


0.05±0.03 


0.05±0.03 







^F(H/3)=(4.112±0.004)xl0-i=^ erg s-^ cm-^. 

^Corrected for extinction with C(H/3)=0.225 and EW(abs)=0.6A. 

'^Corrected for extinction with C(H/3)=0.195 and EW(abs)=0.4A. 
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Table 5. Fluxes of optical and NIR emission lines in Mrk 930 



CLOUDY 



Ion 


100xF(A)/F(H^)^ 


100x/(A)//(H/3)b 


100x/(A)//(H/3)^ 


stars 


stars+shock 




a) Optical spectrum 








0.373 [0 ii] 


207.01ib3.37 


250.9o±4.94 


249.57ib4.97 


259.92 


261.91 


0.384 H9 


3.66±0.31 


9.16ib0.92 


10.30ibl.04 


7.70 


7.68 


0.387 [Ne iii] 


29.83±0.70 


35.40±1.29 


35.30ibl.30 


38.61 


39.17 


0.389 He i+HS 


12.66ib0.46 


20.00ibl.03 


21.19ibl.07 


21.48 


21.40 


0.389 [Ne iii]+H7 


20.54ib0.58 


28.98zbl.l9 


30.20ibl.23 


27.94 


28.12 


0.407 [S ii] 


1.92zb0.46 


2.20ib0.56 


2.21zb0.56 


1.45 


1.63 


0.410 H(5 


19.58±0.55 


26.82±1.12 


27.19=bl.l4 


26.29 


26.24 


0.434 H7 


39.27±0.82 


46.33ibl.42 


46.28ibl.44 


47.32 


47.29 


0.436 [0 III] 


4.84ib0.38 


5.21zb0.51 


5.26ib0.52 


5.71 


5.84 


0.447 He i 


000 1 1^ 

3.38ib0.38 


3.55ib0.46 


3.60ib0.47 


4.14 


4.12 


0.466 [Fe iii] 


1.71±0.35 


1.73±0.38 


1.76±0.39 


1.69 


1.69 


0.469 He 11 


0.83=b0.23 


0.84ib0.25 


0.85zb0.25 


0.00 


0.80 


0.486 H/3 


100.00±1.76 


100.00±1.82 


100.00±1.85 


100.00 


100.00 


0.496 [0 III] 


147.18±2.50 


140.29ib2.75 


143.91zb2.85 


148.71 


148.74 


0.501 [0 III] 


445.01zb7.23 


420.63ib7.18 


431.96±7.47 


447.61 


447.70 


0.588 He I 


13.09ib0.43 


ll.02ib0.59 


ll.49ib0.62 


10.26 


10.21 


0.630 [0 i] 


7.67±0.41 


6.19±0.46 


6.49±0.48 


3.58 


5.18 


0.631 [S III] 


1.64zb0.48 


1.32ib0.40 


1.39zb0.42 


1.61 


1.61 


0.636 [0 i] 


2.12zb0.26 


1.70±0.26 


1.78=b0.27 


1.14 


1.65 


0.656 Ha 


358.66±5.86 


284.llib5.33 


284.18±5.43 


285.23 


285.75 


0.658 [N 11] 


18.19ib0.51 


14.30ib0.66 


15.04ib0.69 


9.51 


9.78 


0.668 He I 


4.02ib0.30 


3.13ib0.32 


3.30ib0.34 


2.90 


2.88 


0.672 [S 11] 


34.74±0.75 


26.97=b0.92 


28.41zb0.97 


9.47 


10.31 


0.673 [S 11] 


25.53±0.65 


19.79±0.78 


20.85±0.83 


6.71 


7.30 


0.707 He I 


2.46ib0.22 


1.85±0.24 


1.95ib0.25 


2.21 


2.20 


0.714 [Ar III] 


9.27ib0.34 


6.92ib0.44 


7.33zb0.47 


6.93 


6.92 




b) Near infrared spectrum 








0.953 [S III] 


61.39ibl.l7 


38.30ibl.l6 


36.77ibl.l2 


38.48 


38.60 


0.955 Pae 


4.68±0.32 


3.77±0.34 


3.82zb0.34 


3.68 


3.67 


1.005 Pa5 


7.78±0.26 


5.81±0.37 


5.84±0.37 


5.50 


5.49 


1.029 [S 11] 


1.71±0.15 


2.00ib0.25 


2.15zb0.27 


0.25 


0.25 



-30- 



Table 5 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H/3)^ 


100x/(A)//(H/3)b 


100x/(A)//(H/3)" 


stars 


stars+shock 


1.083 He I 


36.44±0.65 


22.43±0.80 


21.87±0.78 


20.00 


19.90 


1.094 Pa7 


14.16±0.30 


9.21±0.46 


9.10±0.46 


8.84 


8.82 


1.257 [Fe ii] 


3.42±0.19 


2.23±0.31 


2.19±0.30 


1.30 


1.74 


1.282 P&P 


27.07±0.47 


15.81±0.64 


15.51±0.63 


15.81 


15.79 


1.643 [Fe II] 


3.70±0.21 


1.97±0.20 


1.93±0.20 


1.26 


1.69 


1.736 Brio 


2.14±0.12 


1.45±0.17 


1.50±0.17 


0.93 


0.93 


1.876 Paa 


69.09±3.63 


36.36±2.14 


35.82±2.12 


32.28 


32.25 


1.944 Br 6 


3.70±0.20 


2.14±0.21 


2.15±0.21 


1.79 


1.79 


2.058 He i 


1.73±0.09 


1.09±0.14 


1.12±0.14 


1.48 


1.47 


2.122 H2 1-0 S(l) 


1.13±0.06 


0.76±0.12 


0.79±0.12 






2.165 Br7 


4.59±0.24 


2.50±0.22 


2.50±0.23 


2.68 


2.67 


2.223 H2 1-0 S(l) 


0.85±0.05 


0.46±0.09 


0.43±0.09 






2.248 H2 2-1 S(l) 


0.65±0.04 


0.35±0.08 


0.32±0.08 







^F(H^)=(6.937±0.051)xlO-i^ erg s'^ cm-^. 
^Corrected for extinction with C(H/3)=0.310 and EW(abs) 
'^Corrected for extinction with C(H/3)=0.270 and EW(abs): 



)=2.85A. 
=3.55A. 



Table 6. Fluxes of optical and NIR emission lines in Mrk 996 



CLOUDY 

Ion 100xF(A)/F(H^)'^ 100x/(A)//(H/3)b 100x/(A)//(H/3)^ sta^s 



a) 



0.373 


[0 II] 


149.06±2.33 


0.377 


Hll 


0.83±0.29 


0.380 


HIO 


1.54±0.34 


0.384 


H9 


2.59±0.28 


0.387 


[Ne iiij 


38.04ib0.73 


0.389 


He I+H8 


11.82±0.38 


0.389 


[Ne iii]+H7 


18.47±0.45 


0.403 


He I 


2.24±0.26 


0.407 


[S ii] 


1.12±0.16 


0.410 


Hd 


20.60±0.45 


0.434 


H7 


45.21±0.73 


0.436 


[0 III] 


13.27±0.34 


0.447 


He I 


6.35±0.23 


0.466 


[Fc III] 


2.40±0.19 


0.486 


H/3 


100.00±1.49 


0.492 


He I 


1.07±0.17 


0.496 


[0 III] 


98.73±1.47 


0.501 


[0 III] 


304.25±4.44 


0.520 


[N I] 


0.62±0.11 


0.588 


He I 


19.00±0.32 


0.630 


[0 I] 


4.40±0.12 


0.631 


[S III] 


2.46±0.11 


0.656 


Ha 


377.58±5.49 


0.658 


[Nil] 


26.68±0.41 


0.668 


He I 


5.78±0.14 


0.672 


[S II] 


22.06±0.34 


0.673 


[S II] 


17.44±0.28 


0.707 


He I 


10.71±0.19 


0.714 


[Ar III] 


10.86±0.19 


0.732 


[0 II] 


4.27±0.10 


0.733 


[0 II] 


2.79±0.09 



spectrum 



198.11ib3.78 


199.88±3.80 


157.30 


6.01±2.77 


6.07±2.80 


4.64 


7.14±2.04 


7.21±2.06 


5.95 


8.37±1.21 


8.46±1.22 


7.98 


49.01±1.41 


49.55±1.41 


44.15 


20.60±0.98 


20.84±0.99 


19.70 


28.33±1.08 


28.69±1.08 


29.87 


2.77±0.38 


2.80±0.38 


2.79 


1.37±0.24 


1.39±0.24 


2.16 


29.95±1.07 


29.49±1.06 


26.99 


55.46±1.39 


54.88±1.38 


47.94 


14.83±0.68 


15.13±0.68 


17.46 


6.86±0.45 


7.02±0.46 


6.10 


2.45±0.28 


2.51±0.28 


2.36 


100.00±1.57 


100.00±1.58 


100.00 


1.01±0.20 


1.04±0.20 


1.52 


92.55±1.74 


95.58±1.79 


79.20 


281.75±4.42 


291.21±4.57 


238.40 


0.55±0.13 


0.57±0.13 


0.33 


14.88±0.56 


15.55±0.58 


18.62 


3.24±0.25 


3.40±0.26 


1.62 


1.81±0.19 


1.90±0.19 


3.89 


270.24±4.53 


270.07±4.54 


280.93 


18.91±0.61 


19.89±0.64 


25.85 


4.04±0.27 


4.26±0.28 


4.15 


15.35±0.54 


16.17±0.56 


14.39 


12.11±0.47 


12.76±0.50 


11.22 


7.11±0.36 


7.51±0.37 


13.34 


7.14±0.35 


7.55±0.37 


6.69 


2.74±0.21 


2.90±0.22 


5.36 


1.78±0.17 


1.89±0.18 


4.26 
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Table 6 — Continued 



CLOUDY 



Ion 


100xF(A)/F(H^)^ 


100x/(A)//(H/3)b 


100x/(A)/7(H/3)^ 


stars 


0.775 [Ar iilj 


2.32±0.09 


1.40±0.15 


1.49±0.1o 


1.61 


n nn'7 To txtI 

0.907 [d IIIJ 


23.82±0.40 


13.41±0.47 


1 A /I r» 1 r» c n 

i4.4U±U.oU 


1 A An 

i4.4(J 


n nco rc tttI 
0.95o [d IIIJ 


70.07±1.07 


35.50±0.90 




00.71 




b) Near infrared spectrum 






0.953 [h IIIJ 


70.07±3.91 


35.50±2.12 


38.20±2.29 


35.71 


0.955 rae 


3.85±0.40 


2.56±0.31 


2. 7d±0. 33 


3.67 


1.005 F&S 


9.29±0.54 


5.17±0.36 


5.58±0.38 


5.46 


1.083 He I 


132.86±7.32 


62.56±3.71 


67.64±4.01 


218.83 


1.094 Pa7 


15.32±0.85 


7.70±0.49 


8.33±0.52 


8.82 


1.257 [Fe II] 


2.30±0.15 


l.OlibO.ll 


1.10±0.12 


2.53 


1.279 He I 


2.12±0.18 


0.93±0.11 


1.01±0.12 


0.53 


1.282 Pa/3 


30.25zbl.67 


13.60±0.83 


14.78±0.90 


15.54 


1.643 [Fe lij 


2.28±0.13 


0.92±0.10 


l.OOibO.ll 


2.44 


1.736 Brio 


2.50±0.16 


1.23±0.12 


1.35±0.13 


0.92 


1.944 Br(5^ 


4.82±0.38 


2.07±0.19 


2.27±0.21 


1.74 


1.957 Ha 1-0 S(3) 


0.36±0.03 


0.15±0.05 


0.16±0.05 




2.034 H2 1-0 S(2) 


0.28±0.03 


0.11±0.04 


0.12±0.04 




2.122 H2 1-0 S(l) 


1.06±0.06 


0.40±0.06 


0.44±0.07 




2.165 Br7 


6.07±0.34 


2.46±0.18 


2.70±0.20 


2.60 


2.248 H2 2-1 S(l) 


0.34±0.03 


0.13±0.04 


0.14±0.04 




2.355 H2 2-1 S(0) 


0.21±0.03 


0.08±0.03 


0.08±0.03 





^F(H/3)=(7.241±0.023)xl0-i^ erg s'^ cm-^. 

^Corrected for extinction with C(H/?)=0.450 and EW(abs)=2A. 

'^Corrected for extinction with C(H^)=0.420 and EW(abs)=2A. 



-33- 



Table 7. Fluxes of optical and NIR emission lines in SBS 0335-052E 



CLOUDY 



Ion 


W0xF{X)/F{Rl3)'' 100x/(A)//(H^)'^ 


100x/(A)//(H/3)^ 


stars 


stars+shock 




a) Optical spectrum 








0.373 [0 ii] 


28.34±0.47 


31.11±0.91 


31.43±0.90 


27.36 


27.17 


0.375 H12 


2.25±0.18 


c\ f\ A in on 

2.94±0.39 


2.83±0.39 


3.43 


3.38 


0.377 Hll 


3.85±0.16 


1 n 1 n r' r\ 

5.10±0.60 


4.89±0.60 


4.26 


4.20 


0.380 HIO 


4.99±0.17 


6.15±0.52 


6.01±0.52 


5.52 


5.44 


0.384 H9 


7.61±0.18 


8.98±0.56 


8.88±0.56 


7.49 


7.38 


0.387 [Ne iii] 


22.27±0.36 


c\ A no 1 n ^0 

24.23ib0.78 


24.52±0.78 


24.66 


23.73 


0.389 He I+H8 


16.03=b0.28 


1 T n d 1 n f*c\ 

17.94±0.69 


18.00±0.68 


19.11 


18.80 


0.389 [Ne in]+H7 


23.59±0.37 


25.99±0.82 


26.19±0.82 


23.38 


23.35 


0.403 He I 


1.30±0.13 


1.40±0.20 


1.42±0.20 


1.72 


1.66 


0.407 [S li] 


n n 1 n "1 

0.80±0.13 


0.86±0.16 


0.87±0.16 


0.12 


0.13 


0.410 H5 


24.19±0.37 


c\f* no 1 n 01 

26.28±0.81 


25.78±0.80 


25.81 


25.46 


0.434 H7 


45.88±0.67 


A 0/^ 11 in 

48.3d±1.12 


47.79±1.11 


47.29 


47.09 


0.436 [0 III] 


9.82±0.17 


1 n no in Ad 

10.23±0.48 


10.45±0.48 


7.98 


7.87 


0.447 He i 


3.72±0.12 


3.84±0.29 


3.93±0.30 


3.57 


3.45 


0.469 He 11 


c\ c\r\ 1 r\ r\r\ 

2.20±0.08 


2.22ib0.22 


2.29±0.22 


0.00 


2.21 


0.471 [Ar iv]+He I 


1.29±0.08 


1.30±0.17 


1.34±0.17 


1.00 


1.03 


0.474 [Ar iv] 


1.18±0.07 


1.19±0.16 


1.22±0.16 


0.37 


0.40 


0.486 up 


100.00±1.44 


100.00±1.46 


100.00±1.46 


100.00 


100.00 


0.492 He i 


1.08±0.09 


1.07±0.16 


1.11±0.16 


0.95 


0.92 


0.496 [0 III] 


101.28±1.46 


100.12±1.72 


103.51±1.77 


89.07 


86.16 


0.501 [0 III] 


304.90±4.36 


300.32±4.34 


310.74±4.58 


268.10 


259.34 


0.588 He I 


11.76±0.19 


11.02±0.47 


11.53±0.49 


9.02 


8.70 


0.656 Ha 


300.67±4.29 


273.73±4.38 


273.80±4.38 


282.51 


284.85 


0.658 [N 11] 


0.81±0.04 


0.74±0.12 


0.78±0.13 


0.65 


0.70 


0.668 He i 


3.07±0.06 


2.78±0.23 


2.93±0.24 


2.44 


2.35 


0.672 [S 11] 


2.39±0.06 


2.16±0.21 


2.27±0.21 


0.63 


0.80 


0.673 [S 11] 


2.12±0.05 


1.92±0.19 


2.02±0.20 


0.48 


0.60 


0.707 He I 


5.21±0.09 


4.64±0.30 


4.91±0.31 


3.30 


3.12 


0.714 [Ar III] 


1.93±0.05 


1.71±0.18 


1.82±0.19 


2.07 


1.98 


0.728 He i 


0.79±0.04 


0.70±0.12 


0.74±0.12 


0.69 


0.66 


0.901 PalO 


1.39±0.07 


1.20±0.15 


1.28±0.16 


1.83 


1.80 
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Table 7— Continued 



CLOUDY 



Ion 


100xF(A)/F(H^)'^ 


100x/(A)//(H/3)b 


100xI(A)//(H^)^ 


stars 


stars-hshock 


0.907 [S III] 


3.66±0.09 


3.09±0.24 


3.30±0.25 


3.07 


2.94 


0.953 [S III] 


8.17±0.22 


6.72±0.35 


7.24±0.37 


7.61 


7.29 


0.955 Pae 


3.86±0.19 


3.22±0.25 


3.45±0.27 


3.48 


3.43 




b) Near infrared spectrum 








0.953 [S III] 


8.17±0.58 


6.72ib0.52 


7.24±0.55 


7.61 


7.29 


0.955 Pae 


4.33±0.36 


3.63±0.34 


3.89±0.36 


3.48 


3.43 


1.005 PaS 


5.99±0.41 


4.88±0.39 


5.25±0.41 


5.19 


5.11 


1.083 He I 


32.78±2.14 


26.36±1.78 


28.53±1.93 


27.61 


26.56 


1.091 He I 


0.49±0.06 


0.39±0.08 


0.42±0.09 


0.25 


0.24 


1.094 Pa7 


9.26±0.61 


7.47±0.54 


8.08±0.58 


8.29 


8.17 


1.279 He I 


1.14±0.09 


0.90±0.13 


0.98±0.14 


0.46 


0.44 


1.282 Pa/3 


15.75±1.03 


12.40±0.86 


13.49±0.93 


14.81 


14.72 


1.641 Brl2 


0.63±0.05 


0.50±0.09 


0.54±0.10 


0.54 


0.53 


1.643 [Fe II] 


0.27±0.04 


0.21±0.07 


0.23±0.08 


0.31 


0.82 


1.681 Brll 


0.64±0.05 


0.51±0.09 


0.55±0.10 


0.67 


0.66 


1.736 Brio 


0.97±0.07 


0.75±0.11 


0.82±0.12 


0.86 


0.85 


1.876 Pact 


37.39±2.52 


28.42±2.00 


31.15±2.19 


29.11 


29.05 


1.944 BiS 


2.42±0.17 


1.84±0.19 


2.02±0.20 


1.62 


1.59 


2.058 He i 


1.76±0.12 


1.33±0.15 


1.46±0.17 


1.37 


1.30 


2.122 H2 1-0 S(l) 


0.21±0.02 


0.16±0.05 


0.17±0.06 






2.165 Br7 


3.60±0.24 


2.72±0.24 


2.99±0.26 


2.40 


2.36 


2.248 H2 2-1 S(l) 


0.14±0.04 


0.11±0.08 


0.11±0.08 







''F(H/?)=(6.672±0.012)xl0-i'^ erg s"! cm-^. 

^Corrected for extinction with C(H^) =0.135 and EW(abs)=0.5A. 

^Corrected for extinction with C(H/3)=0.105 and EW(abs)=0.35A. 
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Table 8. Extinction-corrected fluxes of strong hydrogen lines 



Galaxy 



Line 


II Zw 40 


Mrk 71 


Mrk 930 


Mrk 996 SBS 0335- 


-052E 


Theor.^ 








a) Optical spectrum 






0.410 R6 


32.0 


26.8 


26.8 


29.5 


25.8 


26.2 


0.434 H7 


48.9 


48.7 


46.3 


54.9 


47.8 


47.3 


0.486 H/3 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


0.656 Ha 


283.9 


278.7 


284.1 


291.2 


273.8 


279.0 








b) NIR spectrum 






1.094 P7 


7.9 


9.4 


9.2 


8.3 


8.1 


8.6 


1.282 P/3 


13.2 


15.9 


15.8 


14.8 


13.5 


15.2 


2.165 Br7 


2.8 


2.4 


2.5 


2.7 


3.0 


2.5 



^Recombination ratios for rp=15000 K and Np 



Hummer Storey I (|l987l ) 



100 cm" 



from 



Table 9. Extinction coefficients C(H/3) derived from individual hydrogen line to H/3 flux 

ratios^ 



Line 






Galaxy 






II Zw 40 


Mrk 71 


Mrk 930 


Mrk 996 


SBS 0335-052E 








a) Optical spect 


rum 




0.410 m 


0.686±0.215 


0.112±0.150 


0.348±0.287 


0.188±0.224 


0.051±0.159 


0.434 H7 


1.001±0.247 


0.091±0.216 


0.416±0.317 


0.005±0.244 


0.026±0.221 


0.656 Ha 


1.181±0.115 


0.217±0.112 


0.348±0.129 


0.403±0.114 


0.093±0.112 








b) NIR spectrum 




1.094 P7 


1.066±0.180 


0.231±0.180 


0.341±0.069 


0.381±0.184 


0.048±0.217 


1.282 P/3 


1.036±0.164 


0.200±0.163 


0.331±0.052 


0.387±0.166 


0.016±0.197 


2.165 Br7 


1.167±0.139 


0.163±0.139 


0.305±0.136 


0.443±0.141 


0.175±0.168 



^T he adopted theoretical rec ombination ratios are for Te=15000 K and A'^e 
from iHummer &: Storey I (jl987l ) (see Table [8]) . 



100 cm 



-3 
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Table 10. H 


2 emission- 


-line flux ratios 






Line 






Galaxy 


Model^ 


II Zw 40 


Mrk 71 


Mrk 930 


Mrk 996 SBS 0335-052E 


Fluor 


Coll 


1.162 H2 


2-0 S(l) 




O.libO.l 






0.4 


0.0 


1.207 H2 


3-1 S(2) 


0.2±0.1 


0.2±0.1 






0.2 


0.0 


1.233 H2 


3-1 S(l) 


0.4±0.2 


0.4±0.2 






0.5 


0.0 


1.238 H2 


2-0 Q(l) 


0.6±0.3 


0.3±0.2 






0.4 


0.0 


1.311 H2 


4-2 S(l) 


0.4±0.2 


0.3±0.2 






0.4 


0.0 


1.314 H2 


3-1 Q(l) 


0.4±0.2 


0.3±0.2 






0.6 


0.0 


1.601 H2 


6-4 Q(l) 


0.3±0.1 


O.libO.l 






0.4 


0.0 


1.957 H2 


1-0 S(3) 


0.7±0.2 


0.9±0.3 




0.4±0.1 






2.034 H2 


1-0 S(2) 


0.4±0.1 


0.2±0.1 




0.3±0.1 


0.5 


0.3 


2.073 H2 


2-1 S(3) 


0.3±0.1 


0.2±0.1 






0.2 


0.0 


2.122 H2 


1-0 S(l) 


1.0±0.2 


1.0±0.2 


1.0±0.2 


1.0±0.2 1.0±0.4 


1.0 


1.0 


2.223 H2 


1-0 S(0) 


0.4±0.1 


0.5±0.2 


0.6±0.2 




0.6 


0.3 


2.248 H2 


2-1 S(l) 


0.6±0.2 


0.5±0.2 


0.5±0.1 


0.3±0.1 0.6±0.4 


0.5 


0.0 


2.355 H2 


2-1 S(0) 


0.2±0.1 


0.2±0.2 




0.2±0.1 


0.3 


0.0 



'^Model values are from iBlack &: van Dishoeck I (|l987l ) . We adopt their model 1 for fluorescent 
lines and model SI for collisionally excited lines. 
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Table 11. Input parameters for stellar photoionization CLOUDY models 



Param(^(;(\r 


11 Zw 10 


Mrk 71 


Mrk 930 


Mrk 996'' 


SBS 0335-052E 


Q(H) 


53.136 


52.064 


53.893 


52.750 


53.562 




50000 


50000 


50000 


50000 


50000 


Ne, cm ^ 


100 


100 


10 




100 


log / 


-1.10 


-0.30 


-2.30 


-4.58 


-1.60 


log He/H 


-1.07 


-1.07 


-1.05 


-0.93 


-1.09 


log C/H 


-4.68 


-4.68 


-4.30 


-4.47 


-5.57 


log N/H 


-5.29 


-5.83 


-5.34 


-4.90 


-6.13 


log 0/H 


-4.00 


-4.15 


-3.92 


-3.97 


-4.66 


log Ne/H 


-4.71 


-4.87 


-4.70 


-4.67 


-5.40 


log Si/H 


-5.35 


-5.35 


-5.35 


5.57 


-6.05 


log S/H 


-5.63 


-5.78 


-5.70 


-5.50 


-6.53 


log Cl/H 


-7.49 


-7.59 


-7.49 


-7.47 




log Ar/H 


-6.20 


-6.41 


-6.37 


-6.34 


-7.06 


log Fe/H 


-5.80 


-5.90 


-5.99 


-5.70 


-6.10 



^For Mrk 996, a two-zone model for the electron number density spatial 
distribution has been adopted. For the inner high-excitation zone with radii 
r between 0.8 pc and 80 pc, the electron number density has a power-law 
dependence A^e ~ r^"^, decreasing from 4.5x10^ cm~^ to 450 cm~^. For 
the lower-excitation outer zone with r> 80 pc, the electron number density 
has been set to a constant value of 450 cm~^. 
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Table 12. Fluxes of NIR [Fe ll] emission lines relative to Br7 



CLOUDY 

Ion F{X)/F{Bij) /(A)//(Br7)'^ /(A)//(Br7)'^ stars stars+shock 

a) II Zw 40 



1.257 


[Fe 


n] 


0.076 


0.112 


0.109 


0.149 




1.643 


[Fe 


n] 


0.073 


0.088 


0.087 


0.146 














b) Mrk 71 






1.257 


[Fe 


II] 


0.034 


0.036 


0.037 


0.108 


0.159 












c) Mrk 930 






1.257 


[Fe 


II] 


0.745 


0.892 


0.876 


0.485 


0.652 


1.643 


[Fe 


II] 


0.806 


0.788 


0.772 


0.597 


0.633 












d) Mrk 996 






1.257 


[Fe 


II] 


0.378 


0.411 


0.407 


0.973 




1.643 


[Fe 


II] 


0.376 


0.374 


0.370 


0.938 












e) SBS 0335-052E 






1.643 


[Fe 


II] 


0.075 


0.077 


0.077 


0.121 


0.347 



^Corrected for extinction derived from the observed Balmer decrement. 

'^Corrected for extinction derived from the Balmer decrement, itself corrected for 
colhsional excitation of hydrogen Unes. 



- 39 - 



Table 13. Fluxes of MIR emission lines 



CLOUDY 

Ion 100x/(A)//(H/3) stars stars+shock 







a) II Zw 40'^''' 




in K^ 
iU.Oi 


[o IVJ 


40.77zb0.88 


68.19 




1 O C 1 


TAT^ ttI 

[iNe IIJ 


1.31±0.07 


0.43 




io.oo 


FAT^ tttI 

[iNe IIIJ 


24.99±0.66 


33.56 




io. / i 


FG tttI 
[b IIIJ 


ll.25zb0.44 


10.40 




o on 

zo.yy 


\Ti^^ ttI 

[re IIJ 


0.33zb0.09 


0.13 








b) Mrk 930=''^ 




12.81 


[Ne ii] 


5.11zb0.37 


3.14 


3.35 


15.55 


[Ne III] 


18.12zb0.17 


27.82 


28.06 


25.99 


[Fe ii] 


2.70zb0.62 


0.45 


0.65 






c) Mrk 996'''-^ 




10.51 


[S IV] 


5.15zb0.02 


4.02 


4.948 


12.81 


[Ne ii] 


5.76zb0.05 


7.20 


7.26S 


15.55 


[Ne III] 


16.32zb0.76 


15.50 


16.81S 


18.71 


[S hi] 


ll.21zb0.76 


9.98 


10.54s 


25.89 


[0 IV] 


1.82zb0.45 


0.00 


1.428 


33.48 


[S III] 


16.44zb7.27 


12.99 


13.818 






c) SBS 0335-052Eh'^ 




10.51 


[S IV] 


9.22zb0.51 


5.33 


5.27 


15.55 


[Ne III] 


7.99zb0.45 


7.52 


7.28 



^Spitzer data fro m IWu et al. I (j2006l ) and 



O'Halloran et al. I (|2006l b 



^Observed H/3 fl ux of 3.45x10 erg s ^ cm ^ from 



Lagos et aL I ([20071). 



^Spitzer data from lO'Halloran et al. I (|2006l ). 

'^O bserved flux of (2.54zb0.10) x 10"^^ erg s"^ cm^^ 
from iMoustakas Sz Kennicutt I (j2006l ) . 



"Spitzer data from iThuan et al. I (|2008l ). 
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^Extincti on-corrected H/3 flu x of 6.6x10 erg s ^ 
cm-2 from lxhuan et al. I ((2003). 



SModel D from lThuan et al. I (120081 ). 

^Spitzer data from lHouck et al. I (|2004l ). 

'Observed fluxes normalized to the extinction- 
corrected B.0 flux of (1.813±0.G27) x 10"^^ erg s"^ cm"^ 
(jlzotov et al. Il2006bl ). 
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Fig. 1. — 3.5m APO/TripIeSpec NIR spectrum of II Zw 40 in five orders. In each panel, the 
noisy regions of the upper spectrum are omitted. They are caused by insufficient sensitivity 
or strong telluric absorption. The flux scale on the y-axis corresponds to the upper spectrum. 
The lower spectrum is downscaled by a factor of 50 as compared to the upper spectrum. It 
is shown for the whole wavelength range in each order. 
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Fig. 3. — Same as in Fig. [H but for Mrk 930. 
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Mrk 996 
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Fig. 4. — Same as in Fig. [H but for Mrk 996. 




Fig. 5.— 



Wavelength (/lim) 
Same as in Fig. [H but for SBS 0335-052E. 
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Fig. 6. — Redshift-corrected spectral energy distribution of II Zw 40 in the 0.36 /^m - 2.4 
jj,m wavelength range 
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Fig. 7. — Same as in Fig. |6l but for Mrk 71. 
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Mrk 930 
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Fig. 8. — Same as in Fig. |6l but for Mrk 930. 
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SBS 0335-052E " 
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Fig. 10.— Same as in Fig. El but for SBS 0335-052E. 



